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^ ■ ABSTRACT 

b 

, Context. BL Lacertae objects are variable at all energy bands on time scales down to minutes. To construct and interpret their spectral energy 
distribution (SED), simultaneous broad-band observations are mandatory. Up to now, the number of objects studied during such campaigns is very 
limited and biased towards high flux states. 

Aims. We intend to study the simultaneous broad-band emission of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 
1ES 2344+514 by means of a pre-organised multi- wavelength campaign independently of the current flux state. Long-term observations at dif- 
ferent energies might allow for a discussion of the results in a more general context. By correlating the obtained light curves and describing the 
corresponding SEDs with state-of-the-art models we intend to derive physical parameters of the emission region. 

Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were 
complemented by Metsahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term 
F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using 
a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. 

Results. 1ES 2344+514 is found at very low flux states in both X-rays and very high energy gamma rays. Variability is detected in the low fre- 
quency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was potentially caused by shock acceleration 
characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over 
time scales of eleven years, opposite to previous findings . A significant correlation between the 15.0 GHz and the R-band monitoring light curves 
cannot be established. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE en- 
ergies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED 
modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. 
Consequently it remains unclear whether a so-called quiescent state was found in this campaign. 1 
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1. Introduction 

The number of known extragalactic very high energy (VHE, 
^ lOOGeV) gamma-ray sources has been increasing steadily in 
the past seven years and now exceeds 50 (November 2012fl 
Most of these sources are X-ray bright BL Lacertae (BL 
Lac) objects. In BL Lacs the relativistic jet points close to 
our line of sight and the resulting large relativistic beaming 
causes rapid variability in all energy regimes from radio wave- 
lengths to VHE gamma rays. The spectral energy distribu- 
tion (SED) of these objects shows two peaks; the low energy 
peak is attributed to synchrotron emission, emitted by relativis- 
tic electrons spiralling in the magnetic field lines of the jet, 
while the high energy peak is generally considered to be pro- 
duced by inverse Compton scattering. The seed photons for 
the Compton scattering can be the synchrotron photons them- 
selve s (synchrotron self Comp ton, SSC, e.g. Maraschi et alJ 
[1993 M 

oom & Marscherl [l996l) or photons from an external 
radiation fiel d (accretion disk, broad fine region clouds or in- 
frared torus: iDermer & Schlickeiserl 119931: ISikora et al.1 Il994t 
iBIazejowski et al.l2000l) . An alternate source has been proposed, 
that the gamma rays are produced by hadronic processes, that 
is by proton initiat ed cascades or directly through proton syn- 
chrotr on radiation (Man nheim & Biermannl fl992l : iMiicke et al.l 
120031) . 

BL Lac objects were historically divided into two sub- 
classes, depending on the energy of the synchrotron peak. The 
class boundaries can be loosely defined such that low energy 
peaking BL Lac objects (LBLs) have their peak at 10 14 ~ 15 Hz 
(optical regime) and high energy peaking BL Lacs (HBLs) 
at > 10 15 Hz (UV to hard X-rays) (e.g. iPadovani & Giommil 
1995). The class intermediate to these two was introduced by 
Laurent-Muehleisen et alJ d 1999b . noting that BL Lacs exhibit a 
continuous range in SED peak energy rather than a dichotomy. 
The BL Lac sources detected in VHE gamma rays mostly 
belong to the HBL class. Their spectral energy distributions 
can be described with one-zone SSC emission, but the mod- 
elling requires rather high jet speeds while Very Long Baseline 
Interferometry (VLBI) observations have shown that the parsec- 
scale jets of these objects are comparably slow (Lorentz fac - 
tor r mo d e i « 25 compared to Tvlbi ^ 5: iPiner etalJ 12010). 
There fore it has been suggested that the jet is decelerating 
dGeor ganopou los & Kaz anas 2003]) or has a spine-and-sheath 
structure dGhisellini et alJl2005h . Recent VLBI observations of 
the electric vector position angle and fractional polarisation 
distribution in Te V blazars support the spine-sheath scenario 
(IPiner et alJl2010h . 

BL Lac objects show variability at all bands from radio to 
VHE gamma rays. The variability amplitudes vary between the 
different energy regimes and from source to source. The VHE 
gamma-ray detected X-ray selected BL Lacs are typically quite 
faint and mildly variable in the radio, show a large range of 
variability in the optical band and are strongly variable in X- 
rays. In the gamma-ray band they are often mildly variable at 
sub-GeV-GeV energies, while in VHE gamma rays some of 
the sources show extreme variability partly with amplitudes ex- 
ceeding one order of magnitude and flux doubling time scales 
as short as minutes (e.g. Mrk 421, Mrk 501, PKS 2155-304 ; 
lAcciari et al ] l2011at lAlbert et aLlkOPTat lAharonian etai1l2007l) 
whereas others vary with smaller amplitude (e.g. 1ES 1215+303, 
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PG 1553+1 13: lAleksic etaD l2012allbh . The variability is typi- 
cally ^e^cribedjri terms of "quiescence" and "flaring" epochs 
(e.g. lAcciari et al.l 1201 Id) . The quiescent state is characterised 
by a low flux and little variability at all wavebands while the 
flaring epochs exhibit a high flux and large variability. 

Due to their variability and their broad-band emission, the 
SED of BL Lacs has to be based on simultaneous observa- 
tions at all energy ranges (simultaneous multi-wavelength [MW] 
campaigns). For many sources the observations are concen- 
trated on flaring epochs due to a higher detection probability. 
Simultaneous MW observations from radio to VHE gamma rays 
in low flux states were for a long time scarce for these objects 
due to limited sensitivity of the first generation of gamma-ray in- 
struments. Even today such observations are mostly available for 
the three brightest objects, Mrk 42 1 , Mrk 5 01 and PKS2 155 -3 04 
(see e.g. the most rece nt campaigns in Abd o et al.l 1201 1 alibi 
lAbramowski et alJl2012h . 

1ES 2344+5 14 is an HBL at redshift z = 0.044 
dPerlman et al.ll 1996b - It was first detected at VHE gamma rays 
(above 300 GeV) by the Whipple telescope in 1995 during a 
flare with a flux F (> 3 50 GeV) = (6.6 + 1.9) • 10" phcm 2 s 1 
(ICatanese et al.l 119981) and was only the third known extra- 
galactic VHE gamma-ray source. Follow-up observations in 
a lower state did not result in detections with high statisti- 
cal significance unti l the MAGIC observations in 2005-2006 
(lAlbert et al.l l2007bl). The s ource was not seen by EGRET 
(e.g. Mukheriee et al] 1 19971) but was detected by the Fermi 
space telescope with a flux F(l - 100 GeV) = (1.55 ± 0.18) • 
10~ 9 phcirT 2 s _1 and a hard power law spectral index (1.72 ± 
0.08) a s reported in the F erm/-LAT Second Source Catalog 
(2FGL; iNolan et al l 1201 2|) (see also Sect. |4~33T >. Like most 
HBLs it does not exhibit strong variability in the Fermi band 
(variability index ~ 28 in 2FGL, while an index of > 41 was re- 
quired to reject t he null hypothesis of no variability at the 99 % 
confidence level: INolan et al.ll2012l) . Note that 1ES 2344+514 is 
formally not listed as a "clean" source in the Fermi AGN Catalog 
due to its low galactic latitude but nevertheless appears in the 
corresponding source tables. 

In the X-ray band the source is br ight with a 2 ke V flux den- 
sity of 1.14//Jy dPerlman et al.llT996h and showed strong spec- 
tral variability with the synchrotron p eak shifting to higher en- 
ergies when the flux level increases dGiommi et alJ 120001) . In 
the high state, the synchrotron peak frequency was at or above 
lOkeV, making 1ES 2344+5 14 one of the few so-called "ex- 
treme blazars" dCostamante et al.ll200ll) with synchrotron peak 
frequencies in the hard X-rays. Chandra observations revealed 
diffuse X-ray emiss ion as well as seven individual point sources 
in its environment ( Dona to et al . 2003}. 

In the optical band the overall brightness of the source shows 
only very moderate variability (of the order of 0.1 mag). This is 
due to the bri ght host galaxy wh ich contributes ~ 90 % to the 
observed flux dNilsson et al]l2007l) . 

In the radio band the source is rather faint with a core 
flux density S core (5 GHz) aa 0.07 Jy measured by VLBI 
dGirolettietal] 12004) and an overall flux density on arcsecond 
scales of Sarcsec (5 GHz) = (0.23 + 0.01) Jy (average of 18 F- 
GAMMA0 single-dish observations from 02/2007 to 04/2009). 
Using Very Long Baseline Array (VLBA) imaging the apparent 
jet speeds of different components have been determined to be 
$ 3 c with the most robust me asurement of (0.62 ± . 05) c found 
for o ne individual feature dPiner & Edwards! 12004b IPiner et aD 
2010). The lower frequency Very Large Array (VLA) maps 
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Fig. 1. Sky map of the region around 1ES 2344+514 (marked 
by a green cross). Green radio contours at 1.4 GHz are over- 
laid on an R-band image. IR point sources are indicated by red 
circles. An "X" and its label mark indi vidual components iden- 
tified by Chandra dDonato et al.l l2003h . The logarithmic grey 
scale is showing scaled densities. Radio contours are given from 
0.001 Jy/beam until 0.241 Jy/beam in 20 logarithmically scaled 
steps. Only IR sour ces of J magnitude < 15 are displayed. Data 
reference: X-rays: iDonato et all d2003l) : optical: DSS2red; IR: 
2MASS; radio: NVSS, obtained from NED. 



(kpc scale) showed an extended and complex radio structure at 
1 .4 GHz with ~ 45° misalign ment compared to higher frequency 
(5 GH z, pc scale) radio maps dRector et al.ll2003b iGiroletti et al.l 
12004 . 

The combination of archival, non-simultaneous data in the 
radio, optical and X-ray regime reveals that only one of the indi- 
vidual X-ray components in the field of view of 1ES 2344+514 
is bright at 1.4 GHz (component "E", see Fig. [TJ). This com- 
ponent coincides v ery well with the radio feature reported by 
iRector et all (120031) and IGiroletti et al.l (12004 . but is not present 
in the IR or R-band. Consequently, there are no other potential 
VHE candidate sources in the immediate vicinity of the source 
at an angular separation smaller than the MAGIC angular reso- 
lution of ~ 0.1°. The nature of the radio feature can not be iden- 
tified unambiguously. Pulsars, being faint i n the optical regime , 
would still be viable candidates. However. IGiroletti et al ] (12004 
found a connection of the emission between the feature and the 
core in VLA radio images. Also the proximity between these two 
(angular distance of ~ 180", i.e. ~ 160 kpc) indicates that they 
might be related. The jet of the AGN may bend on kpc scales 
by ~45° and interact with the intergalactic medium, resulting in 
a radio hot spot. The wide opening angle of the jet and the low 
surface brightness on these scales do not support the interpreta- 
tion of the feature as a hot spot at this distance from the core 
though. Moreover, this would be in contradiction to the unifica- 
tion scenario whe re the BL Lacs are sugge sted to be beamed FR- 
I radio galaxies (IUrrv&Padovani| ["l995). Note, howeve r, tha t 
similar resul t s have been found by e.g. lLandt & Bignali (|2008); 
iKharb et al l d2010h . Future VLBI measurements of the radio 
spectrum of the feature may distinguish between the radio hot 
spot or foreground/background source interpretation. 



Table 1. Multi-wavelength observations of 1ES 2344+514. The 
dates are given in MJD-54700 and rounded down. In the case of 
OVRO, RATAN-600, KVA+Tuorla, Swift XRT and MAGIC-I, 
the given observation periods were not covered continuously. 



Instrument 


Band" 


Observation Date 


Effelsberg 


radio 


56; 78; 106; 155 


IRAM 


radio 


25; 106; 137; 178 


Metsahovi 


radio 


30; 46; 124; 127; 130; 136; 138 


OVRO 


radio 


61-179 


RATAN-600 


radio 


29-42 


VLBA 


radio 


61-62 


CrAO 


R-band 


74; 77; 85; 101; 105; 112; 117 


KVA+Tuorla 


R-band 


22-134 


Swift 


UV & X-rays 


30; 45-84 


AGILE 


HE gamma rays 


70-100 


Fermi 


HE gamma rays 


59-100 


MAGIC-I 


VHE gamma rays 


59-100 



To date, 1ES 2344+514 has been studied in only one MW 
campaign that included gam ma-ray observations , conducted by 
RXTE Swift and VERITAS dAcciari et al.l201 lbh . lGiommi et al.l 
d2012l) reported about Planck, Swift and Fermi observations, 
covering energies from radio to GeV, but detecting the source 
only in the UV and X-ray bands. In this paper we present the 
first simultaneous radio to VHE gamma-ray observations of 
1ES 2344+514. The campaign was organised independently of 
the flux state to make investigations of a low, possibly the "quies- 
cent", state of the source possible. The observations were sched- 
uled to give the best simultaneous coverage between the different 
instruments, with less than a day time difference between VHE, 
X-ray and optical bands. The time delays with respect to radio 
observations were larger due to the longer variability time scale 
in this energy regime. The campaign took place in autumn 2008 
shortly after the launch of the Fermi gamma-ray satellite. In to- 
tal six radio observatories contributed, including VLBA imag- 
ing of the source in several frequency bands. 1ES 2344+5 14 was 
monitored in the optical R-band by the CrAO, KVA and Tuorla 
telescopes, in ultraviolet and X-rays by Swift UVOT and XRT 
and in high energy (HE) gamma rays by AGILE and Fermi. The 
core part of the campaign was centred around the MAGIC VHE 
gamma-ray observations of the source. Parts of t he MW data sets 
have been presented in Rtig amer et al.l d201 1 alibi) . In this paper 
we present the complete results of the MW campaign. 

The paper is organised as follows: in Sect. 2, we present short 
descriptions of the various participating instruments, their obser- 
vations as well as the corresponding data analyses. The results 
will be shown in Sect. 3 and discussed in Sect. 4 including the 
spectral energy distributions and the theoretical models. Final 
remarks are given in Sect. 5. 



2. Instruments, Multi-Wavelength Observations and 
Data Analysis 

In this section, the instruments participating in the MW cam- 
paign, their observations and data reduction processes will be 
presented ordered by their wavelength regime. A summary of 
the observation dates is given in Table [TJ 



Notes. <a) The exact energy bands are given in Sect. [2] 
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2.1. The MAGIC Telescope 

The MAGIC (Major Atmospheric Gamma-ray Imaging 
Cherenkov) project operates a system of two 17-m Imaging Air 
Cherenkov Telescopes loca ted on the Canary Isla nd of La Palma 
2200 m above sea level (lAleksic et al.l 1201 2d) . MAGIC has 
been operating in stereoscopic mode since 2009, accordingly 
the observations presented in this paper were conducted with 
MAGIC-I only (mono mode). MAGIC-I had a standard trigger 
threshold of 60GeV for observations at low zenith angles, an 
angular resolution of ~ 0.1° for single events and a n energy res- 
olutio n above 150 GeV of ~25% (for details, see lAlbert et al.l 
120081) . 

MAGIC-I observed 1ES 2344+514 from 20/10/2008 to 
30/11/2008 at zenith angles between 23° and 31° for a total 
of 26.4 hours in so-called wobble mode, where the source was 
displaced by 0.4° from the camera center in order to allow the 
recording of simultaneo us OFF-source dat a with the same offset 
from the camera center (Daumetal 1 119971) . 

The data were analysed as described in lAleksic et al.l d2010l) 
with the exception of the signal arrival time extraction. Instead 
of determining the arrival time of the signal at the pulse maxi- 
mum, which was needed at that time due to the special nature 
of those data, the standard method of determining the signal ar- 
rival time at half of the rising flank was used here. 20.8 hours of 
data survived the quality selection. Background suppression was 
accomplished by a cut in shower area versus shower SIZE (i.e. 
total photoelectron content), optimised on 0.7 hours of data from 
a high state of Mrk421 taken during the same observing period 
as 1ES 2344+514 and hence with similar data quality and obser- 
vation conditions. The significance of the signal was determined 
by a cut in (f optimised also on the Mrk421 data set, where 6 
is the angular distance between the expected and reconstructed 
source position. All significances of the VHE signals given in 
the fo llowing sections were determined by Eq. 17 of iLi & Ma 
(1983) with a — 1/3, i.e. using 3 OFF regions. 

The source spectrum has been derived from events with 9 2 < 
0.046 deg 2 , yielding an analysis threshold of ~ 190 GeV. Upper 
limits (UL) were calculated by applying model 4 of iRolke et alJ 
(120051) using a confidence level (c.l.) of 95 %. The conversion 
from the differential spectrum to spectral energy density vF v has 
been accomplished by multiplying the differential flux with the 
energy of the Lafferty-Wyatt bin center (lLaffertv & Wvatl l995) 
squared. 

The MAGIC analysis results presented here were confirmed 
by an independent internal analysis. 

2.2. The AGILE Satellite 

AGIL E (Astrorive latore Gamma ad Immagini LEggero) 
dTavani et al.l 12009) is a scientific mission of the Italian Space 
Agency dedicated to the observation of astrophysical sources of 
high energy gamma rays in the 30 Me V - 50 GeV energy range, 
with simultaneous X-ray imaging capability in the 18-60keV 
band. AGILE is the first high-energy mission which makes use 
of a silicon detector for the gamma ray to pair conversion. The 
AGILE payload combines for the first time two coaxial instru- 
ments: the Gamma-Ray Imaging Detector (GRID, composed 
of a 12-planes Silicon-Tungsten tracker, a Cesium-Iodide mini- 
calorimeter and the anti-coincidence shield) and the hard X-ray 
detector Super-AGILE. The use of the silicon technology pro- 
vides good performance of the gamma ray GRID imager in a 
relatively small and compact instrument: an effective area of the 
order of 500 cm 2 at several hundred MeV, an angular resolution 



of around 3.5° at 100 MeV, decreasing below 1° above 1 GeV, 
a very large field of view (~2.5sr) as well as accurate timing, 
positional and attitude information. 

During the period 07/2007- 10/2009, AGILE was operated 
in "pointing observing mode", characterised by long observa- 
tions called Observation Blocks (OBs), typically of two to four 
weeks duration, mostly concentrated along the Galactic plane. 
Since 1 1/2009 the satellite has been operating in "spinning ob- 
serving mode", surveying a large fraction (about 70 %) of the sky 
each day. The time period covered by the 2008 MW campaign 
includes the AGILE OB 6400, publicly available from the AS DC 
Multimission Archive web page0. 1ES 2344+514 was observed 
by AGILE at ~ 40° off-axis from the mean pointing direction in 
the time window 31/10/2008 to 30/1 1/2008. 

AGILE-GRID data from the official Processing 
Archive (SPINNING sw = 5_21_18_19 and POINTING 
sw = 5_19_18_17), obtained by u sing the AGILE Standard 
Analysis Pipeline (Pitto ri et al.l 120091) . were analysed using the 
latest scientific software (AGILE_SW_5.0_SourceCode) and 
in-flight calibrations (10023) publicly available since 30/09/201 1 
at the ASDC siufl Counts, exposure, and galactic background 
gamma-ray maps were created with a bin-size of 0.3° x 0.3°, 
for E > 100 MeV, selecting only events flagged as confirmed 
gamma-ray events. Events collected during passages of the 
South Atlantic Anomaly or whose reconstructed directions form 
angles with the satellite-Earth vector smaller than 90° were 
rejected to avoid Earth albedo contamination. In order to derive 
the estimated flux (or flux upper limits) of the source we ran the 
AGILE point source analysis software based on the maximum 
likelihood technique using a radius of 10°. 

2.3. Fermi-LAT 

The Fermi sat ellite started with the official science data tak- 
ing 4/08/2008 dAtwood et al.H 2009). Two different detectors are 
on board: the Gamma-ray Burst Monitor (GBM), sensitive at 
low energies (8 keV-40MeV), and the Large Area Telescope 
(LAT), sensitive at high energies (20MeV-> 300 GeV). 

Typically, the Fermi satellite is rocked first towards the north 
pole of the orbit and then, in the next orbit, towards the south, 
alternating in this way the pointing in every orbit. This main 
operating mode, called "All-Sky scanning mode", allows for full 
sky coverage every two orbits, or three hours. 

The LAT is a large field of view (~2.4sr) electron-positron 
pair conversion telescope made up of a high-resolution silicon 
microstrip tracker, a Csl hodoscopic electromagnetic calorimeter 
and an anti-coincidence detector for the identification of charged 
particle backgrounds. The full d escription of the instr ument and 
its performance can be found in lAtwood et al] (120091) . The LAT 
point spread function (PSF) depends strongly upon the energy 
of the impinging gamma ray and on the depth of the conversion 
point in the tracker, and mildly upon the incidence angle. For 
1 GeV normal-incidence conversions in the upper section of the 
tracker the PSF 68 % containment radius is 0.8°. 

The Fermi-hPX data for 1ES 2344+5 14 presented here were 
obtained in the time period between 20/10/2008 22:35:00 UTC 
and 30/11/2008 21:31:00 UTC coordinated with the observa- 
tions with MAGIC. The data have been analysed by using the 
standard Fermi-LPtl Science Tools software package, version 
09-27-01 as described in the Cicerone websitfl The Pass 7 

3 http://www.asdc. asi.it/mmia/index.php ?mission=agilemmia 

4 http://agile.asdc.asi.it/publicsoftware.html 

5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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Source e vent class and P7SO URCE_V6 instrument response 
functions (lAtwood et al.l2009l) were used in our analysis. We se- 
lected events in a region of interest (Rol) centred on the source 
position within 15°, having an energy between 100 MeV and 
300 Ge V. In order to avoid background contamination from the 
bright Earth limb, time intervals when the Earth entered the LAT 
Rol were excluded from the data set. In addition, events with 
zenith angle s larger than 100° with respect to the Earth refer- 
ence frame (Abdo et al. 2009a) were excluded from the analy- 
sis. The data were analysed wi th an unbinned max imum like- 
lihood technique, described in iMattox et al.l ([1996), using the 
analysis software (gtlike) developed by the LAT team and de- 
scribed in the Cicerone website mentioned above. The fitting 
procedure maximises the likelihood acting simultaneously on 
the free spectral parameters for the source of interest, those of 
nearby gamma-ray sources and the diffuse backgrounds, mod- 
elled using ringJ2year_P76ji'0 for the Galactic diffuse emission 
and isotropJ2year_P76jourcejv0 for the extragalactic isotropic 
emission modelfl To maintain comparability, photon fluxes 
were converted to spectral energy densities applying the same 
method as used for AGILE. 

In addition we also performed a dedicated analysis of 
the highest energy photons (> 100 GeV) detected from 
1ES 2344+514 within the first 44 months of LAT operation. 
Only events of the purest class (Pass_7_V6_Ultraclean) from a 
68 % containment radius around the direction of the source were 
considered for this anal ysis. Front and bac k photons, accord- 
ingly to the definition in Atw ood et al.l d2009l) . were treated sepa- 
rately, having a different distribution of the PSF. Since no results 
on such events over this long time scale have been reported in 
literature, the analysis has been applied to four additional TeV 
HBLs with a comparable redshift (Mrk421, Mrk501, Mrk 180 
and 1ES 1959+650). 

2.4. Swift 

The Swift satellite ( iGehrels et al.ll200"4T) is equ ipped with three 
telesc opes, the Burst Alert Telescope (BAT; IBarthelmv et al.l 
120051) which covers the 14-195keV en ergy range, the X- 
ray telescope (XRT; iBurrows et a II 12001 covering the 0.2- 
lOkeV energy band , and the UV/Optical Telescope (UVOT; 
iRoming et alj 120051) covering the 180-600nm wavelength 
range with V, B, U, UVW1, UVM2 and UVW2 filters. 

Swift XRT observed 1ES 2344+514 from 09-1 1/2008 with 
a total of 21 exposures (see Table lA. Il l with exposure times rang- 
ing from 200 s to 5 ks. The two exposures lasting well below 1 ks 
were too short for deriving a flux and were therefore excluded 
from the analysis. The XRT data were processed with stan- 
dard procedures using the FTOOLS task XRTPIPELINE (ver- 
sion 0.12.6) distributed by HEASARC within the HEASOFT 
package (v.6.10). E vents with grades 0-12 were selected (see 
IBurrows et al.l2 005) and latest response matrices available in the 
Swift CALDB (v.20 100802) were used. For the spectral analy- 
sis the source events were extracted in the 0.3-10keV range 
within a circle with a radius of 20 pixels (~47"). The back- 
ground was extracted from an off-source circular region with a 
radius of 40 pixels. The spectra were extracted from the cor- 
responding event files and binned using GRPPHA to ensure a 
minimum of 25 co unts per energy bin, in order to guarantee re- 
liable x 2 statistics dGehrelsl fl986). Spectral analyses were per- 
formed using XSPEC version 12.6.0. The spectral index was de- 
termined using an absorbed power law fit {fyE~ a -e~ T ) from 0.3 - 



lOkeV, with the absorption r being the product of the absorp- 
tion hydrogen-equivalent column density Nh and the element- 
specific energy-dependent photoelectric cross section cr (E). Nh 
was fixed to the G alactic value in the d irection of the source 
of 1.5 ■ 10 21 cm" 2 dKalberlaet al.ll2005l) . The analysis yielded 



a value of (2.0 + 0.2) • 10 cm . Since some daily data sets 
showed hints of spectral curvature, also fits using a log-parabola 
model (/o • E~( a+bl ° Sl >> ( - E ^ . e ~ T ) were performed. However, the 
log-parabola fit was significantly preferred by a logarithmic like- 
lihood ratio test over the simple power law model (see Table 
Table |A~TT > in only a minority of the fits. Therefore, as a com- 
mon basis the simple power law results were used. 

For the long-term source evolution, 67 observations of 
1ES 2344+514 between 2005 and 2010 were analysed. A 
slightly different analysis procedure was used compared to 
the MW data reduction. The spectra were determined using 
XSELECT (V2.4b) to extract events with an energy of 0.5- 
lOkeV from the corresponding event files. The background was 
deduced from an annulus around the source with an inner radius 
of 50 pixels (~ 118") and an outer radius of 70 pixels (~ 165"). 
Spectral analysis and binning was performed in ISIS (V 1.6.2- 
3), where a minimum signal to noise ratio of 5 was required 
for grouping the data. The spectral index was determined in the 
range 0.5- lOkeV using an absorbed power law fit. To calcu- 
late the integral flux the photon flux was evaluated on a fine 
grid between 2 and 1 ke V. The neutral hydrogen-equivalent col- 
umn density was determined for each spectrum from the spec- 
tral fit, yielding for spectra with a d.o.f. > 35 a mean value of 
(1.71 ± 0.14) • 10 21 crrT 2 . Flux errors are given at a 90 % confi- 
dence level. The event counts for calculating the hardness ratios 
for the MW data were extracted applying this pipeline in the full 
energy range. 

Swift UVOT observed the source with all filters (V, B, U, 
UVW1, UVM2, UVW2). The source counts were extracted from 
a circular region 5 arcsec-sized centred on the source posi- 
tion, while the background was extracted from a larger circu- 
lar nearby source-free region. These data were processed with 
the uvotmaghist task of the HEASOFT package. The ob- 
served magnitudes have been corrected for Gala ctic extinction 
Eb-v = 0.191 mag (Schlaflv & Finkb einerll201 lh using the ex- 
tinction curve from iFitzpatrickl d 1999b adopting Ry = 3.07 
dMcCall & ArmouriEz OOO)'. The host-galaxy flux contributes sig- 
nificantly to the observed brightness in the V-, B- and U-bands, 
however no values for the contribution were found in the lit- 
erature. Therefor e, the contrib u tion i s estimated from the R- 
band value from iNilsson et all ( 20071) ("apertu r e 5") using the 
galaxy colors at z = from iFukugita et al.l d 1995b resulting 
in V = (1.96 + 0.16) mJy, B = (0.95 ± 0.20) mJy and U = 
(0.22 + 0.20) mJy. In these bands the host galaxy contributes 
~ 80 - 90 % to the measured flux and additionally the uncertainty 
of the host-galaxy contribution is rather large. Therefore these 
bands are not considered for spectral energy distribution mod- 
elling. 

The magnitudes measured in the UV filters were converted to 
un its of ergcm~~ 2 s " 1 usi ng the photometric zero points as given 
in iBreeveld et al.l d201 ll) and effective wavelengths and count- 
rate-to-flux ra tios of GRBs from the Swift UVOT CALDB 02 
(v.20101 130V lRaiteri et alJ d2010l) noted that these ratios are not 
necessarily eligible for BL Lac objects, due to their different 
spectrum and a B-V value typically larger than the applica- 
ble range. Therefore, they determined the UVOT effective wave- 
lengths and count-rate-to-flux ratios anew (for BL Lacertae, an 



LBL at z = 0.069). We compare these values with the ones used 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html in this work and find that the difference amounts to < 1 % for 



5 



J. Aleksic et al.: The Simultaneous Low State SED of 1ES 2344+514 from Radio to VHE 



the V, B and U filters. In the case of the UV bands, the effective 
wavelengths (count-rate-to-flux ratios) are ~ 7 % (~ 2 %), ~ 3 % 
(~ 1 %) and ~9 % (~ 13 %) larger for the UVW1, UVM2 and 
UVW2 filters, respectively. Therefore we did not apply a new 
calibration but conservatively increased the error of the UVW2 
count-to-flux ratio from ~ 2 % to 13 %. However the actual un- 
certainty should be much below that, considering that some (if 
not most) of the difference between the ratios arises solely from 
using new effective wavelengths, which is not the case in this 
work. 

Swift BAT operates in full sky mode. The BAT data of 
1ES 2344+514, taken from the 58-Month Catalog^, have been 
re-binned using the tool rebingausslc from the HEASOFT 
package to weekly (7 days), monthly (30.44days), quarterly 
(91.31 days) and yearly (365.24days) bins. The default settings 
for the bin centre of rebingausslc have been used, no trials 
have thus been made for sel ecting the binning. Integral fluxes 
were calculated according to Tuel ler et alj j2010) by multiplying 
the Crab-normalised count rate of 1ES 2344+514 with the Crab 
flux measured in the same time interval and energy band. These 
fluxes were then converted to spectral energy densities in each 
energ y band at the Lafferty-Wy att bin position ( Laffertv & Wvatt 
1 19951) assuming a simple power law with a spectral index of 2.62 
as given in the BAT 58-Month Catalog. 

2.5. KVA and Tuorla 

1ES 2344+514 has been monitored in the optical R-band by the 
Tuorla Blazar Monitoring Program since 20020 The observa- 
tions are done using the Tuorla 1-m telescope (Finland) and 
the Kungliga Vetenskapsakademien (KVA) 35-cm telescope (La 
Palma). The latter can be controlled remotely from the Tuorla 
Observatory. In the following, "KVA" will be used as a syn- 
onym for "KVA+Tuorla". The source is typically observed a 
few times per week, but during the Swift pointings mechanical 
problems prevented KVA observations. The photometric mea- 
surements are made in differential mode, i.e. by obtaining CCD 
images of th e target and calibrate d comparison stars in the same 
field of view dFiorucci et al.ll9 98). The magnitudes of the source 
and comparison stars are measured using aperture photometry 
and the (colour corrected) zero point of the image determined 
from the comparison star magnitudes. The object magnitude is 
computed using the zero point and a filter-dependent color cor- 
rection. Magnitudes are then transferred to linear flux densities 
using the formula F = Fo ■ 10 mag/ ~ 25 , where mag is the magni- 
tude of the object and Fo is a filter-depende nt zero point (in the 
R-band the value Fo = 3080 Jy is used from lBessellll 19791) . 

Since 1ES 2344+514 has a bright host galaxy and a nearby 
star that contributes to the observed flux, these contributions 
have to be removed in o rder to derive the co re flux for the spec- 
tral energy distribution. Nils son et alJ (120071) determined these 
contributions which depend on seeing and the aperture used for 
the measurement. Since all observations for this campaign were 
done with constant aperture (7.5") and in similar seeing condi- 
tions, we subtract a constant value of (3.70 + 0.05) mJy. 

2.6. CrAO 

Observations from the Crimean Astrophysical Observatory 
(CrAO) were obtained with the AZT- 1 1 telescope and an FLI 

7 After an update, the 58-Month Catalog contains as of now 
(10/2012) the results from the first 66 months of observation. 

8 http://users.utu.fi/kani/lm/ 



IMG 100 IE CCD camera, through an R-band filter. Differential 
photometry was performed between the blazar and published 
comparison stars on the same CCD frame. The comparison stars 
and apertures used were the same as for KVA. The resulting 
magnitudes were converted to mJy using the standard formula. 
The CrAO flux densities were found to be ~ 12% lower than 
the KVA points and were shifted by a fixed value (~0.49mJy) 
to match the KVA observations. The corresponding shift has 
been deduced from the average flux density difference between 
both telescopes for nights with an observation time difference 
< 0.3 days. Two out of seven data point pairs satisfied this con- 
dition. A difference of ~ 10 % is expected due to CrAO using the 
Johnson R-band filter whereas KVA is measuring in the Cousins 
R-band filter. 

2.7. Effelsberg 100-m and IRAM 30-m Radio Telescopes 

Quasi-simultaneous cm-to-mm radio spectra have been ob- 
tained within the framework of a Fermi related monitoring 
progra m of gamma-ray blazars, namely the F-GA MMA pro- 
gram dFuhrmann etalj|2007l lAngelakis etal] I2008I) . The total 
frequency range spans from 2.64GHz to 228.4GHz using the 
Effelsberg 100-m and IRAM 30-m telescopes. The millimetre 
observations are closely coordinated with the more general flux 
monitoring conducted by IRAM, and observations of both pro- 
grams are included in this paper. 1ES 2344+514 has been ob- 
served in autumn 2008 once a month with these facilities. 

The Effelsberg measurements were conducted with the sec- 
ondary focus heterodyne receivers at 2.64, 4.85, 8.35, 10.45, 
14.60, 23.05, 32.00 and 43.00GHz. The observations were per- 
formed quasi-simultaneously with "cross-scans" (that is, slew- 
ing over the source position in azimuth and elevation direction), 
with an adaptive number of sub-scans for reaching the desired 
sensit ivity (for details see lFuhrmann et alJ20 08: Angelakis et al.l 
2008). Consequently, pointing off-set corrections, gain correc- 
tions and atmospheric opacity corrections have been applied to 
the data. The conversion to Jy has been done using the standard 
calibrators: 3C48, 3C 161, 3C286, 3C295 and NGC7027. The 
standard deviation of the flux calibrators amounts to < 5 % at 
43.00 GHz and < 1 % at 2.64 GHz. The Effelsberg error bars are 
given including systematic uncertainties. 

IRAM (Institut de Radioastronomie Millimetrique) operates 
a 30-m radio telescope located on Pico Veleta near Granada 
in Spain. The IRAM observations of 1ES 2344+514 and pri- 
mary/secondary calibrators were carried out with calibrated 
cross-scans using the receivers operating at 86.2 and 142.3 GHz, 
occasionally also at 228.4GHz. The opacity corrected scans 
were converted into the standard temperature scale and finally 
corrected for small remaining pointing offsets and systematic 
gain-elevation effects. The conversion to the Jy flux density scale 
was done using the instantaneous conversion factors derived 
from the frequently observed primary (Mars, Uranus) and sec- 
ondary (W3(OH), K3-50A, NGC 7027) calibrators. 

2.8. Metsahovi 14-m Radio Telescope 

The 37 GHz observations were conducted with the 13.7-m di- 
ameter Metsahovi radio telescope, which is a radome-enclosed 
paraboloid antenna in Finland. The measurements were made 
with a 1 GHz-band dual beam receiver centred at 36.8 GHz. The 
HEMPT (high electron mobility pseudomorphic transistor) front 
end operates at room temperature. The observations are ON- 
ON observations, alternating the source and the sky in each feed 
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horn. A typical integration time to obtain one flux density data 
point is between 1200 and 1400 s. The detection limit of the tele- 
scope at 37 GHz is of the order of 0.2 Jy under optimal condi- 
tions. Data points with a signal to noise ratio < 4 are considered 
as non-detections. 

The flux density scale is set by observations of DR 21. The 
sources NGC 7027, 3C 274 and 3C 84 are used as secondary cal- 
ibrators. A detailed description of the data reduction and analysis 
is given in Ter asranta et al.l (Q998). The error estimate in the flux 
density includes the contribution from the measurement rms and 
the uncertainty of the absolute calibration. 



2.9. OVRO 40-m Radio Telescope 

Regular 15.0 GHz observations of 1ES 2344+514 were carried 
out as part of a high-cadence gamma-ray blazar monitoring pro- 
gram using th e Owens Valley Rad io Observatory (OVRO) 40- 
m telescope dRichards et alJl201lt) . This program, which com- 
menced in late 2007, now includes about 1600 sources, each ob- 
served with a nominal twice per week cadence. Data during the 
beginning of this MW campaign were unavailable due to a hard- 
ware outage. The OVRO 40-m results used in this paper span the 
period 22/10/2008 to 1 1/02/2012. 

The OVRO 40-m uses off-axis dual-beam optics and a cryo- 
genic high electron mobility transistor (HEMT) low-noise am- 
plifier with a 15.0 GHz centre frequency and 3 GHz bandwidth. 
The total system noise temperature is about 52 K, including re- 
ceiver, atmosphere, ground, and CMB contributions. The two 
sky beams are Dicke-switched using the off-source beam as a 
reference, and the source is alternated between the two beams 
in an ON - ON fashion to remove atmospheric and ground con- 
tamination. A noise level of approximately 3-4mJy in quadra- 
ture with about 2 % additional uncertainty, mostly due to point- 
ing errors, is achieved in a 70 s integration period. Calibration 
is achieved using a temperature-stable diode noise source to re- 
move receiver gain drifts. The flux density scale is derived from 
observations of 3C 286 assuming a value of 3.44 Jy at 15.0GHz 
(Baa rs et al.ll 1977b . The systematic uncertainty of about 5 % in 
the flux density scale is not included in the error bars. Complete 
details of the r eduction and calibration procedure are found in 
iRichardset ai1d2011l) . 



2.10. RATAN-600 

The radio spectrum of 1ES 2344+514 w as observed with the 
600-m ring radio telescope RATAN-600 dKorolkov & Pariiskii 
1 19791) of the Special Astrophysical Observatory, Russian 
Academy of Sciences, located in Zelenchukskaya, Russia, from 
20/09/2008 to 03/10/2008. The continuum spectrum was mea- 
sured six times quasi-simultaneously (within several minutes) 
in a transit mode with six different receivers at the follow- 
ing central frequencies (and frequency bandwidths): 0.95 GHz 
(0.03 GHz), 2.3 GHz (0.25 GHz), 4.8 GHz (0.6 GHz), 7.7 GHz 
( 1 .0 GHz), 1 1 .2 GHz ( 1 .4 GHz), 2 1 .7 GHz (2.5 GHz). Due to ra- 
dio frequency interference, we were unable to detect the source 
at the two longest wavelengths. An average spectrum of the six 
independent 5-22 GHz measurements is presented in this paper. 
Details on the method of observatio n, data process i ng, an d am- 
plitude calibration are described by iKovalev et al.l ([1999). The 
data were collected using the southern sector with the Flat re- 
flector. 



2.11. VLB A 

1ES 2344+514 was observed with the VLBA dNapierll 19951) on 
23/10/2008 at 4.6, 5.0, 8.1, 8.4, 15.4, 23.8 and 43.2 GHz in 
the framework of a surve y of parsec-scale radio sp ectra of 20 
gamma-ray bright blazars dSokolovskv et alj|2010bl) . The obser- 
vations were conducted with ten on-source scans (each four to 
seven minutes long depending on frequency) spread over eleven 
hours. The data reduction w as performed in the standard man- 
ner using the AIPS package (Greisen 1990). A special p rocedure 
similar to the one described in Sokolovsky et al. (201 1) was ap- 
plied to improve amplitude calibration of the correlated flux den- 
sity resulting in ~ 5 % calibration accuracy at the 4.6- 15.4 GHz 
range and ~ 10% accurac y at 23.8 and 43.2 GHz. The Difmap 
software (Shepherd 1997) was used for imaging and modeling 
of the visibility (uv) data. The integrated par sec-scale flux densi- 
ties were derived by summing all CLEAN dHogboml 19741) com- 
ponents used to represent calibrated visibilities. 

1ES 2344+514 was also observed with VLBA at 15.4 GHz 
during the campaign as a part of the MOJAVE (Monitoring 
Of Jets in Active galactic nuclei with VLBA Experiments].] 
long-term program to monitor radio brightness and polarisation 
variations in jets associated with active galaxies visible in the 
northern s ky. The data were analysed using the standard proce- 
dures (see lLister et a l. 2009a b). Elliptical Gaussian components 
were used to determine positions and flux densities of individ- 
ual emission regions within the source. We adopt a cosmology 
with Q. m = 0.27, Q A = 0.73 and H {) = 71 km s 1 Mpc" 1 . 
The MOJAVE archive contains two sets of VLBA data on this 
source at 15.4 GHz. On e set contains four epochs published in 
iPiner & Edwards! d2004l) that span the range 10/1999 to 03/2000. 
The second consists of ten epochs covering 05/2008 to 1 1/2010. 

3. Results 

3.1. Very High Energy Gamma Rays 

The MAGIC data analysis yielded a marginal signal of 3.5 <r 
for the complete data set (see Fig. [2] and for detailed results 
Table IA.2l i. which is below the 5 <x standard for source dis- 
coveries in VHE astronomy. Since 1ES 2344+514 is a well- 
established VHE emitter and the direct environment is lacking 
suitable alternative source candidates (see Sect[T|), we assume 
that the entire excess comes from the source. The rather long 
observation time of ~ 20 hours and the fairly large events statis- 
tics not dominated by individual features in time makes us con- 
fident about the reliability of the signal. Therefore, we derived 
an average spectrum. 

The measured (EBL de-absorbed) spectra are rather well fit- 
ted Or 2 /d.o.f. = 0.36/1 for both of them; see also the residuals 
shown in Fig. [3} by a simple power law of the form 

^ =fo ■ 10- 12 TeV-' cnrV • (E/E y a (1) 
aE 

yielding f = 4.0 + 1.2 (4.8 ± 1.5) at E = 0.5 TeV and a = 
2.4 ± 0.4 (2.2 + 0.4) (see Fig. [3}. The given errors are statistical 
only. We adopt the MAGIC standard systematic errors of 16 % 
on the energy scale, 1 1% on the flux no rmalisation and ±0.2 
on the spectral index dAlbert et al.ll2008l) . The low redshift of 
the source renders differences between the current extragalactic 
background light (EBL) models negligible. Here, the effects of 
EBL absorptio n were corrected by applying the Kneiske "lower 
limit" model dKneiske & Doldl2010l) . 

9 http://www.physics.purdue.edu/MOJAVE/ 
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Significance 3.5o, off-scale 0.33 
217.7 excess events, 2775.3 background events 




AGILE window 
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Fig. 2. 9 2 plot for 1ES 2344+5 14. The vertical dashed line gives 
the 2 cut, defining the signal region. The background events are 
shown by the grey area. 
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Fig. 3. Top panel: Measured (blue, filled circles) and de- 
absorbed (red, open circles) MAGIC spectra for 1ES 2344+514, 
shown togethe r with the MAG IC 2005 spectrum of 
1ES 2344+514 (Albert et al. 2007b) as well as the MAGIC 



Crab spectrum ( Alb ert et al . 2008). Upper limits were derived 
with a c.l. of 95 %. Bottom panel: Fit residuals defined as the 
difference between flux and fit divided by the flux value. 
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Fig. 4. The red points show the daily MAGIC light curve of 
1ES 2344+5 14 derived from this campaign. For points having 
a flux consistent with or below zero, shown by the grey trian- 
gles, 95 % c.l. upper limits were calculated. The red bar on top 
of each upper limit arrow visualises the bin width. A fit with a 
constant to all daily flux points is shown by the red dashed line. 
Additionally, the fluxes for the two observation periods (see also 
Table lA.2l are given as blue full squares. 



3.2. High Energy Gamma Rays 

AGILE did not detect the source during the campaign. The 
AGILE maximum likelihood analysis using the latest in-flight 
calibrations yielded a 95 % c.l. UL on the flux above 100 MeV 
of 3.7 • 10~ 8 phcirT 2 s _1 from an effective exposure of ~ 2.8 ■ 
10 8 cm 2 s for the MW observation period. Searching for short 
flares on time scales of seven as well as two days did not yield 
any detection. Given the non-detection of the source we adopt a 
"standard" spectral photon index of 2. 1 for the likelihood analy- 
sis. 

Fermi-LAT did not detect 1ES 2344+514 between 0.1 and 
300 GeV during the campaign (effective exposure: ~ 3.7 • 
10 9 cm 2 s). The data were searched for short-time variability on 
daily and weekly time scales without a clear sign of such. The 
2-year Catalog public light curve does not show significant vari- 
ability on time scales of months around the time of the MW cam- 
paign. Upper limits at a 95 % c.l. have been determined applying 
the standard Bayesian approach for the MW time slot, assuming 
a spectral index of 2.1 to be consistent with the AGILE calcula- 
tions. These amount to (in ph cm" 2 s" 1 ) 3 .0 ■ 1 0~ 8 (0. 1 - 0. 3 GeV), 
6.7 ■ 10~ 9 (0.3-1.0GeV), 2.7 • 10~ 9 (1.0-3.0GeV), 8.8 ■ lO" 10 
(3.0- lOGeV) and 8.6 • 10~ 10 (10- 100 GeV). 



No significant variability could be found over the entire ob- 
servation period on daily time scales, as can be seen from the 
light curve in Fig.|4](note that the fluxes are calculated subtract- 
ing OFF data from ON data and can therefore become nega- 
tive) and the flux values given in Table IA.2I The overall flux 
F(> 170 GeV) amounted to (7.4 + 2. 1) ■ 10~ 12 ph cm -2 s _1 . A fit 
with a constant yields a ^ 2 /d.o.f. = 19.0/13, which gives ~ 12 % 
probability for a constant flux. The low probability arises domi- 
nantly from the negative fluctuation around MJD 54767 and the 
highest flux point at MJD 54787. The latter is indicating a higher 
state of the source, but since the point is less than 2 <x above 
the fit line, it statistically does not give evidence for variability. 
Fitting the fluxes for the two MAGIC observation periods, the 
^ 2 /d.o.f. = 4.0/1 (probability ~5%), which is still consistent 
with the hypothesis of a constant flux of the source. 



3.3. X-Rays 

Swift XRT detected significant variability (see Fig. [5] and 
Table lA.ll . The 2 - 10 keV flux increased by ~ 50 % within two 
days, followed by a slow decline nearly halving the flux dur- 
ing eight days. In the following the flux rose again, showing 
irregular patterns, and eventually reached the highest flux dur- 
ing these observations on the last day. The quicklook Swift XRT 
intra-day light curves (from the Swift Monitoring Prograrro did 
not show a conclusive hint of intra-day variability during the 
MW campaign. The spectral index derived from a simple power 
law fit between 0.3-10keV is consistent with being constant 
during this campaign (jf 2 /d.o.f. = 10.8/18). On the other hand, 
the evolution of the hardness ratio (defined as the ratio of event 
counts between 2- lOkeV and 0.2-1 keV), another measure of 

10 see http://www.swift.psu.edu/monitoring/ 
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Fig. 5. The combined 1ES 2344+5 14 multi-wavelength light curve derived from this campaign. The observation window of AGILE 
is shown by the vertical solid green lines. Horizontal short-dashed lines represent fits with a constant to the data points. In the 
R-band, the fit is being applied to the combined KVA and CrAO data points. Hardness ratio is defined as the event counts between 
2 - 10 keV divided by the counts between 0.2 - 1 keV. The Swift BAT data are shown for several bin widths some of which have been 
multiplied by a factor for clarity as indicated in the panel legend. The fit shown in the BAT panel has been calculated for the 1 week 
binning. One point from the daily BAT light curve (MJD 54867.5, rate - (6.08 + 3.49) • 10~ 3 cnts s _1 detector -1 ) is not shown for 
clarity of the plot. The red solid line in the MAGIC panel gives the overall measured flux during the campaign. See text for details. 
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the spectral shape, cannot be described satisfactorily by a con- 
stant fit Of 2 /d.o.f. = 31.6/18, see Fig. 0. 

Having found a rather hard spectral index down to ~ 1 .8 in 
the XRT band, the BAT 66 months data have been searched for 
hints of a signal. Indeed, during the time of XRT observations 
there are indications of a positive flux for several consecutive 
days, though insignificant due to limited statistics. Therefore, 
the daily BAT light curve was re-binned to different time scales 
(see Sect. 12.41 ). As can be seen from Fig. |5J variability may be 
present in the weekly binned data (^ 2 /d.o.f. = 30.3/22 for a 
constant flux). The weekly high flux point during the XRT mea- 
surements at MJD 54772.5 has a significance of 4.9. However, 
further analysis shows that this anomalous high flux can be at- 
tributed to an artifact of the BAT coded-mask imaging and hence 
is not believed to be due to any real increase in the emission of 
1ES 2344+514. For the monthly binned points, the probability 
of variability increases (constant flux fit: ^/d.o.f. = 12.8/4). 
The quarterly and yearly results will be discussed in Sect. l4.2l in 
the context of the long-term behaviour of the source. 

3.4. UV and Optical 

KVA found the source on a modest overall flux density level 
of ~4.2mJy when compared to earlier and later KVA measure- 
ments (see also Fig.[T4li. No significant variability throughout the 
entire observation period is found in the R-band. The data points 
are consistent with a constant flux density (^ 2 /d.o.f. = 13.1/29). 
The CrAO points are noisier than the KVA points, but all of them 
are compatible with the KVA data within less than two error 
bars. Applying a constant fit to the combined KVA+CrAO mea- 
surements does not provide evidence for variability (^ 2 /d.o.f. = 
35.8/36). Swift UVOT also did not find significant variability at 
any of the measured frequencies (see Fig.[5]and Table lA.3l . 

3.5. Radio Bands 

3.5.1. Single-Dish Observations 

Single-dish radio observations were conducted from 2.64 GHz 
(Effelsberg) to 228.39 GHz (IRAM). The latter instrument did 
not find a significant signal from the source, hence 3 <x ULs were 
calculated. The measurements conducted by Effelsberg show 
significant variability (although of small amplitude) throughout 
the observations, as can be seen in the top panel of Fig. [5] The 
flux density was rising first towards MJD 54779.0 at all frequen- 
cies (observations at 2.64 GHz were not conducted that day) and 
declined slowly afterwards. RATAN-600 found the source prior 
to the Effelsberg observations on a flux density level consistent 
with the first Effelsberg measurements. The OVRO light curve 
shows no clear evidence for variability, having a probability of 
8.7 % Cr 2 /d.o.f. = 33.9/24) for a constant flux density. When 
omitting the outlier around MJD 54870, the probability for a 
constant flux density is rising to 13.6 %. 1ES 2344+514 was too 
faint to be detected by Metsahovi during the campaign and for 
07/10/2008 (MJD 54746), an upper limit on the flux density at 
37 GHz of < 0.33 Jy with S/N > 4 was calculated. The source 
was detected by Metsahovi three months earlier at a flux density 
level of (0.38 + 0.09) Jy, which is consistent with the derived up- 
per limit. 

In Fig. [6] the quasi-simultaneous (separated by ~five 
days) spectra of Effelsberg+OVRO and VLBA are shown. 
Clearly the VLBA integrated spectrum is much flatter than the 
Effelsberg+OVRO spectrum and can be well fit by a simple 
power law of the form S = v~ ff , where S is the flux den- 



sity. The resulting spectral index a is 0.10 + 0.04. On the con- 
trary, a simple power law (a = 0.42 + 0.01) can not describe 
the Effelsberg+OVRO spectrum sufficientlyOj, judging from the 
residuals in Fig. [6] A broken power law is clearly preferred, 
whose fit applied to the Effelsberg+OVRO data results in the fol- 
lowing parameters: E bKak = (5.6 + 1.0) GHz, a\ = 0.49 + 0.03, 
a 2 = 0.34 + 0.05 and a normalisation of (0.153 + 0.004)Jy at 
10 GHz. 

To understand the radio behaviour of AGNs, they have to be 
studied over long periods of time due to the rather large vari- 
ability time scales compared to e.g. X-rays. 1ES 2344+514 has 
been observed in the past on a regular basis at radio frequencies. 
The combined quasi-simultaneous (time difference < 14 days) 
radio spectra from Effelsberg, Metsahovi, OVRO and RATAN- 
600 from 2007 through 2009 are shown in Fig. [7] (for a time- 
resolved version see Fig. IA. lb . IRAM ULs, where the lowest 
flux density UL is 0.96 Jy, are not shown for clarity. At frequen- 
cies below ~ 20 GHz the source shows steep radio spectra while 
above this frequency, the spectra become flat or inverted. The 
Metsahovi points reveal a high flux density state around MJD 
55039 at 37 GHz with a doubling time of < 28 days and a de- 
cline to the original flux density value of < 15 days. This flare 
was not seen by IRAM two days later (UL at 86.24 GHz of 
1.74Jy, at 142.33 GHz of 1.95 Jy). The quasi-contemporaneous 
OVRO points, from the flaring day as well as 4, 12 or 15 days 
after the flare, did not show a significantly higher flux density. 
However, the flare may have been missed due to the comparably 
sparse sampling during these days. A fit with a constant from 
MJD 55024 until MJD 55054 does not yield significant variabil- 
ity Or 2 /d.o.f. = 4.9/8). There are very few ex amples of such fas t 
variability at 37 GHz for HBLs, e.g. Mrk421 dLichti et al. 2008). 
That is mainly due to their faintness and consequently low de- 
tection rate at this frequency. Nevertheless some of these objects 
are detected at clearly higher flux density values in between pe- 
riods of non-detectio ns, giving a hint for fast variability (see e.g. 
iNieppola et alll2007l) . 

Figure [8] shows the light curve measured by Effelsberg in 
the context of the F-GAMMA program from beginning of 2007 
until mid of 2009 (MJD 54155-54952). Apart from an overall 
higher flux density state especially at high frequencies from the 
start of the observations until mid of 2008 (~ MJD 54600), sev- 
eral structures are visible. These features clearly do not appear 
at all frequencies at the same time but are missing (partly due 
to sampling) or possibly delayed in different bands. The most 
prominent feature is a peak around MJD 54547, which is visible 
at all frequencies where coverage is available. 

3.5.2. Interferometric Observations 

The VLBA image of 1ES 2344+514 (Fig. gj reveals a core- 
dominated structure with a smooth jet extending in South-East 
direction. At the distance of 1ES 2344+514 (z = 0.044), the lin- 
ear scale of the images is 0.9 pc/mas. The integrated parsec-scale 
radio spectrum is flat, which is typical for blazars. The VLBA 
data can be used to estimate the size of the radio core (the un- 
resolved structure at the North- Western end of the jet in Fig. [9]) 
at each frequency (see Table |2]i by modelling it with a circular 
Gaussian emission component. At high frequencies the core size 
can only be constrained to < 0.05 mas, i.e. R < 10 17 cm, while at 
lower frequencies the core is partially resolved. The limiting res- 



11 Note that all error bars shown in Fig.[6]contain the systematic con- 
tribution, because of which x 1 goodness of fits cannot be given. 
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Fig. 6. Top panel: Radio spectrum of 1ES 2344+514 measured 
quasi-simultaneously by Effelsberg, OVRO and VLBA. The 
VLBA points represent integrated flux densities. The solid lines 
illustrate fits with a simple power law, the dashed line shows a 
broken power law fit. Bottom panel: Fit residuals, which are the 
differences between the measured flux densities and the fit val- 
ues divided by the flux density values. 

Table 2. Core size as a function of frequency v. 



V 


Size 


Resolution Limit 


[GHz] 


[mas] 


[mas] 


43.2 


<0.06 


0.06 


23.8 


<0.04 


0.04 


15.4 


0.07 ± 0.04 


0.04 


8.4 


0.14 + 0.05 


0.05 


8.1 


0.09 ± 0.05 


0.06 


5.0 


0.16 + 0.07 


0.08 


4.6 


< 0.11 


0.11 



oluti on and compon e nt size uncertain ties were estimated follow - 
ing lFomalonj (I 1 9991) . iLobanovl (120051) and lKovalev e"t~ai1(l2005h . 

Using multi-epoch MOJAVE results (see Table IA.41 . the av- 
erage core brightness temperature at 15.4 GHz can be deter- 
mined as 7b ~ 8 • 10 10 K. While being rather smooth, the jet 
of 1ES 2344+514 can still be divided into several individual 
emission components. Consistency of their positions, fluxes and 
sizes among MOJAVE epochs suggests that these components 
are real stable structures in the jet, not an artefact of represen- 
tation of a smooth continuous jet with discrete Gaussian com- 
ponents. Analysis of the 15.4 GHz MOJAVE monitoring shows 
no significant motion of the jet components over the entire ob- 
serving period of eleven years. Even across the large eight-year 
time gap, the positional changes of the fitted component posi- 
tions are smaller than their overall scatter in the post-2008 pe- 
riod. Parameters of the jet components and results of the formal 
linear fits to their trajectories are presented in Table lA.4"l Among 
the jet components, C 3 is the brightest and smallest one, located 
~ 0.6 mas from the core. C3 provides the strongest limits on 
the apparent jet speed v app of (-4 ± 7) pas yr~' corresponding 
to y6 app = Vapp/c = -(0.01 +0.02). It can be clearly id e ntified 
with the component C 3 described in iPiner & Edwar ds (2004) 
and IPiner et al] d2010h . where the jS app values derived for C 3 
were given as -0.19 ± 0.40 and 0.10 ± 0.02, respectively. 
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Fig. 7. Radio spectra of 1ES 2344+514 taken by Effelsberg, 
Metsahovi, OVRO and RATAN-600 from 2007 through 2009. 
Data from different instruments have been combined if the time 
difference was less than 14 days. The thin solid lines simply con- 
nect the data points. The legend contains the MJD range of the 
combined spectra. We recall that detection limit of Metsahovi at 
37 GHz under optimal conditions is ~ 0.2 Jy. See also Fig. IA. II 
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Fig. 8. Light curve of 1ES 2344+514 measured by Effelsberg. 
The solid lines connect the points for illustrative purposes. The 
MW campaign duration is visualised by the grey box. 



4. Discussion 

4.1. Individual Energy Bands 

4.1.1. Very High Energy Gamma Rays 



Above 200 GeV, the integral flux amounted to (5.5 + 1.7) • 
10 _12 ph cm 2 s 1 , more than a factor 4 lower than the for- 
mer MAGIC detection (which, at the time, constituted the 
lowest flux measured of this source at VHE). Compared 
to the ave rage flux measured by VERITAS > 300 GeV in 
2008 (see lAcciari et al.l 1201 lbl) . the average flux found here 
(F(> 300 GeV) = (3.4 + 1.0)-10- 12 phcnr 2 s _1 )is still lower by 
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Fig. 9. Naturally weighted CLEAN image of 1ES 2344+5 14 ob- 
tained with the VLBA on 23/10/2008 (MJD 54762) at 15.4 GHz. 
Green circles mark positions and sizes of model components 
listed in the Table lA.4l from North-West to South-East: the Core, 
C3, C2, CI. 1 mas corresponds to 0.9 pc in absolute length. 
Contour image parameters: peak = 0.11 Jy/beam, beam = 0.98 x 
0.58 mas at PA = -42.2°, first contour = 0.50mJy/beam, con- 
tour level increase factor = 2. 



a factor of > 3 and hence represents the lowest flux reported from 
1ES 2344+514 at low VHE thresholds up to now. At high ener- 
gies the HEGRA collaboration reported a flux F (> 970 GeV) = 
(6.0+1.9) • 10 _13 phcm _2 s _1 after 72.5 hours of ob servation 
time between 1997 and 2002 (lAharonian et al.l [20041) . which 
is comparable to the one reported here (F(>970GeV) = 
(4.8 + 3.1)- lO^phcirrV). 

Significant flux variability, as detected by VERITAS down to 
time scales of about one day, is not found within this MW cam- 
paign. Also the MAGIC measurements on monthly time scales 
are statistically consistent with constant, though with a low prob- 
ability. The flux points from individual days are all compatible 
with the mean flux within less than 2 <x. This may on the one 
hand be due to the low flux level and rather large statistical er- 
rors, nevertheless the measurements exclude on the other hand a 
rise in flux by more than a factor of ~ 9 of the mean flux (derived 
from the highest 3 <x UL calculated for all light curve points). 
Note that the integral flux during the VERITAS flare was a fac- 
tor of ~ 20 higher than the one measured in this campaign. 

A spectrum of 1ES 2344+5 14 at VHE has already been mea- 
sured by several instru ments. Form er detections by Wh ipple 
(ISchro edter et all 120051) MAG IC (lAlbert et al.l l2007bl) and 
VERITAS dAcciari et al.ll2.01 lw confirm spectral variability, as 
expected for a BL Lac type object. The spectral index has 
ranged from 2.43 + 0.22 staf + 0.15 SVS f (lAcciari et al.ll201 lbl) to 
2.95 + 0.12 stat + 0.2 syst ^Albert et all |2007bl) with a trend of 
a hardening of the spectrum with increasing flux. Contrary 
to that, the value of 2.4 + 0.4 found here indicates a hard 
spectrum despite a very low flux state. However, these re- 
sults are still consistent with most of the archival measure- 



ments due to the large statistical errors. A hard spectral in- 
dex would imply that the second SED peak was located at 
unusually high energies for that flux level, opposite to the 
spectral hardening trend observed for the best studied blazars 
(e.g. Mrk421, Mrk501, PKS 2155-304: iFossati et al.l [2008; 



Albert et al. 2 0073 lAbramowski et aljfe OlO) 



4.1.2. High Energy Gamma Rays 

AGILE-GRID did not detect the source neither during the MW 
campaign nor over the entire period from 07/2007 up to 01/201 1 . 
For the latter, the AGILE maximum likelihood analysis yielded 
a95%c.l. ULonthe flux > 100MeVof2.7- lO^phcnr 2 s" 1 , 
consistent with the 2FGL average flux above 100 Me V which is 
about 0.9 • lO^phcm^s- 1 . 

1ES 2344+514 is rather dim for a TeV AGN in the Fermi 
band. It wa s detected for the first time after 5.5 months of ob- 
servations dAbdo et al.l l2009bb . LAT did not detect the source 
durin g the MW ca mpaign. From the first (1FGL; Abd o et al.l 
1201 Oh to the second dNolan et al.l2012h LAT Source Catalog list- 
ing, the measured fluxes from 1 - 100 GeV and spectral power 
law indices changed from (1.40 + 0.30) • 10~ 9 phcirT 2 s _1 to 
(1.55 + 0.18) -10- 9 ph cm" 2 s- 1 and 1.57 + 0.12 to 1.72 + 0.08, 
respectively. These values are consistent within the statistical er- 
rors, indicating that the spectral shape did not change signifi- 
cantly on these time scales. Also the monthly light curve shows 
mostly upper limits and marginal detections without signs of ma- 
jor flares. In fact, only one flux point from the monthly binned 
Fermi-LAZX data is available for 1ES 2344+514 within the first 
nine months of regular measurements, the remaining observa- 
tions resulted in ULs. 

Consequently, 1ES 2344+514 seems to be, within the limits 
of the AGILE-GRID and FermZ-LAT sensitivities, a rather sta- 
ble and weak source in the HE gamma-ray band over long time 
scales. Hence, archival data should yield a fairly good estimate 
of the actual flux during this MW campaign. We therefore use 
the spectral information from 1FGL on a quasi-simultaneous ba- 
sis for SED modelling (see Sect. 14. 3t . 

The LAT high energy analysis revealed nine events with en- 
ergies in excess of 100 GeV within the first 44 months of op- 
eration from 1ES 2344+514, the highest energy photon having 
an energy of nearly 500 GeV (see Table We compare these 
with the number of events detected from four similar sources 
(see Sect. I2.31 > in Table I A. 51 An investigation of the distribution 
of event energies is strongly limited by the small event statis- 
tics, but judging from Fig. IA.2I most of them are clustered for 
Mrk421 at 100 GeV, whereas the distribution seems to be shifted 
to ~ 150 GeV forMrk501 and -200 GeV for 1ES 2344+5 14. If 
real, distinct HE flares may be responsible for most of the events 
> 100 GeV detected from Mrk501 and 1ES 2344+514 (we note 
that the events are not clustered in time), in contrast to Mrk421 
for which the distribution indicates a constantly high flux at HE. 

The number of events should be correlated directly with the 
source luminosity. Determining the latter at 60 Ge V (from their 
respe ctive photon fluxes between 10 and 100 GeV in Nola n et all 
2012) and normalising the photon counts to the distance of 
1ES 2344+514, a linear fit for the five sources yields the ex- 
pected correlation with a slope of (0.99 + 0.24) counts per 
10 43 ergs -1 (not shown). This indicates that the 2FGL fluxes 
are a suitable representation of the average source behaviour. 
The goodness of the linear fit is rather low though, having a 
^ 2 /d.o.f. = 7.4/3, but is preferred by a logarithmic likelihood 
ratio test with 98.9% over a fit with a constant (^- 2 /d.o.f. = 
25.5/4). The reason for this are the comparably low number 
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Table 3. Fermi-LNT detected events with an energy > 100 GeV 
within the first 44 months of operation from the direction of 
1ES 2344+514 (R.A. 356.77°, Dec. 51.71°). 



MJD 


Energy 


R.A." 


Dec." 


Sep/ 




[GeV] 


n 


n 


[arcmin] 


54879.961 


221 


356.59 


51.80 


9 


54992.961 


174 


356.79 


51.61 


6 


55041.439 


283 


356.73 


51.75 


3 


55358.826 


495 


356.86 


51.63 


6 


55553.247 


201 


357.01 


51.60 


10 


55896.009 


114 


356.79 


51.68 


2 


55702.733 


207 


356.97 


51.63 


9 


55936.262 


107 


356.79 


51.91 


11 


55948.736 


231 


356.74 


51.73 


2 



Notes. {a) Right ascension (J2000) and {b) declination (J2000) of 
the event. (c) Angular separation between the event direction and 
1ES 2344+514. 
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Fig. 10. Swift XRT spectral index (determined from a simple 
power law fit between 0.3 and lOkeV) versus the integral flux 
for 1ES 2344+514 for this campaign. The dashed line shows a 
fit with a constant, whereas the solid line denotes a linear fit. 



of counts from 1ES 1959+650, which may arise from the flat- 
ter spectral index at HE, and the high number of events detected 
from 1ES 2344+514 (which should be 2-3 according to its lu- 
minosity). Considering the similar luminosities of the sources, 
this should be attributed rather to a higher flaring duty cycle than 
a higher long-term average flux of the source, which would be in 
line with the interpretation of the observed shift in event energy 
distributions. Alternatively, the event counts may also be arti- 
ficially increased by false identification of galactic foreground 
events of 1ES 2344+514, being located at a low galactic lati- 
tude of -9.9°. However, applying the same analysis to two re- 
gions containing no HE source at the same galactic latitude as 
1ES 2344+514, but 2.5° away from the object, did not result in 
the detection of any event with energy > 100 GeV. 

The weakness of this investigation is the low statistical ba- 
sis of only five sources. Additionally, we note that the events 
above 100 GeV have been extracted from 44 months of observa- 
tions, whereas the luminosities were determined from 2FGL (24 
months). These arguments render our conclusions rather specu- 
lative. A catalogue of sources with events > 100 GeV based on 
longer observation times is needed to conduct a more reliable 
study. 

4.1.3. X-rays 

Compared to previous observations, also the soft X-ray flux was 
detecte d at very low levels during this campaign. In lAcciari et al.l 
(I201 lbl) . the lowest reported X-ray fluxes from 2-10keV by 
Swift XRT and RXTE PCA were (9.6 + 0.6) • 10~ 12 erg cm 2 s 1 
and (9.5+2.6) ■ 10~ I2 ergcirT 2 s _1 , respectively. The lowest 
flux in our sample, which was also used to derive the "low 
state" SED (see Sect. 14.3b . is more than 15 % below that level 



((7.9 + 0.5) • 10~ 12 ergctrr 2 s" 1 , see Table|AJ}. Other historical 
measurements show that the source is rather often found in such 
low flux states between 2 and lOkeV (e.g. 8.4- 10~ 12 ergcirT 2 s _1 
measure d by BeppoSAX in 1998, 9 ■ 10~ 12 ergcm" 2 s" 1 by Swift 
in 2005: lGiommi et aLEOOOt ITramacere et al.ll20O7l) . but not be- 
low the lowest flux reported here. 

The XRT spectral index (determined from a power law fit 
from 0.3-10keV) measured during this campaign varied be- 
tween 1.76 ± 0.13 and 2.16 + 0.12, a smaller dynamical range 
compared to previous observations at s imilar energies (see e.g. 
Giommi et al. l2000l lAcriari¥afll2011bh . The index is not sig- 
nificantly variable over time, nevertheless there seems to be a 



dependence on the integral flux, see Fig. [TO] which is generated 
mainly by the highest flux point though. A linear fit results in a 
slope of - (2.77 + 1.1 1)- 10 2 per 10~ 12 ergcuT 2 s _1 with a good- 
ness of fit of 99.9% Ot" 2 /d.o.f. = 4.2/17), whereas a fit with a 
constant has a^ 2 /d.o.f. = 10.8/18 (90.5%). A logarithmic like- 
lihood ratio test prefers the linear fit with 97.9 %. More mean- 
ingful in terms of theoretical models would be to investigate a 
correlation between the spectral index and the peak position, but 
because the latter cannot be determined due to lack of significant 
curvature in the spectra, the integral flux was used. A negative 
correlation between flux and spectral index is expected e.g. for 
an increase of the maximum electron energy in SSC models (e.g. 
iMastichiadis & Kirklll997l) . 

The variable hardness ratio allows to test independently if 
the spectral shape changed considerably with the flux during the 
observations. Especially during high flux states, the hardness ra- 
tio seemed to increase (judging from Fig. |5), meaning that the 
flux at higher energies rose stronger than at lower ones. A weak 
correlation (^- 2 /d.o.f. = 12.6/17 for a linear dependence) be- 
tween the flux and the hardness ratio is visible (see Fig. ITTb . A 
constant fit yields a^ 2 /d.o.f. of 31.6/18. Therefore, according 
to the logarithmic likelihood ratio test, the linear fit is preferred 
with a confidence of 98.9%. This is independently confirming 
the correlation in the index versus flux plane in Fig. [10] and can 
be interpreted as the common "harder spectrum when brighter" 
trend during a blazar flare (see e.g. iPian et aD 19981). which 
has already been reported a lso for 1ES 2344+514 (Giommi et al.l 
2000; lAcciariet al.ll2011bl) . 

Using only the data during the flare, i.e. from MJD 54757 
to 54769, a hint for a counter-clockwise evolution s eems to be 
appare nt in the hardness ratio -flux plane (Fig. [T2l . Kirk et al. 
(Il998l) explain such a behaviour in a model where the flare arises 
from a shock front accelerating electrons within a relativistic jet. 
A counter-clockwise evolution is visible when the observations 
happen close to the maximum emission frequency of the elec- 
trons, where the acceleration and cooling time scales are com- 
parable. In this case, the electrons will not be accelerated to the 
highest energies and no related flare at gamma-ray energies is 
expected. This is in agreement with our simultaneous gamma- 
ray observations, although our VHE light curve does not exclude 
the presence of a flare of similar amplitude to have appeared at 
gamma rays at high confidence. An additional hint for similar 
cooling and acceleration time scales being responsible for the 
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Fig. 11. Swift XRT hardness ratio versus integral flux for 
1ES 2344+514 for this campaign. The dashed line shows a fit 
with a constant, whereas the solid line denotes a linear fit. 



0.7 





























































































































/ 














































































I I I 





i , , 


.ii, 










7 


i 9 1 


1 


1 1 


2 13 14 1 


5 16 



Fig. 12. Swift XRT hardness ratio versus integral flux for 
1ES 2344+514 for the time of the flare (MJD 54757-54769). 
The blue arrows give the evolution in time. 



flare is given by the constant spectral index (x 2 /d.o.f. of 1.4/10) 
during the high flux measurements. 

4.1 .4. Radio Bands 

The inconsistency between contemporaneous spectra of 
Effelsberg+OVRO and VLBA (see Fig. [6j can be explained by 
the different observation techniques. VLBA, due to the inter- 
ferometric measurements, is not sensitive to the steep spectrum 
extended emission from the large scale jet but observes directly 
the flat spectrum of the parsec-scale structure, whereas the 
single-dish telescopes Effelsberg and OVRO measure the whole 
jet. As the brightness of the extended components decreases 
at higher frequencies, the parsec-scale spectrum becomes 
prominent and the single-dish spectrum becomes flatter with 
increasing frequency. The power law indices resulting from 
spectral fitting are in accordance with that interpretation. A 
value of ~ 0.5 is expected in the case of an optically thin kpc 
scale jet, hence a\ ~ 0.49 as found for the Effelsberg+OVRO 
spectrum indicates emission dominated by the large scale jet. 
The flattening or upturn of the spectra at frequencies above 
~ 20 GHz (as seen in Fig. [7]), on the other hand, is not an 
observational artifact but a consequence of high amplitude 
variability of the mm radio emission, originating from a more 
compact region than the one dominating the cm-band radio 
spectrum. 



Fo llowing the classification scheme of Angelakis et al.l 
(I2012I) . the long-term radio spectra shown in Fig. [7] categorise 
1ES 2344+514 as a type 4b object. Their spectra are well de- 
scribed by a simple two component system of a power-law quies- 
cent spectrum (attributed to the optically thin diffuse emission of 
a large scale jet or even a relic of a recent flare or blends of such) 
and a convex synchrotron self-absorbed spectrum (attributed to 
a recent outburst superimposed on the quiescent part). Such an 
outburst may be explain ed in the framework of the model of 
lMarscher&GeaJdl985h where the emission is coming from a 
shock propagating in an adiabatic relativistic jet. The spectra 
from the time of the core of the campaign (MJDs ~ 54757 and 
~ 54779) unfortunately consist only of low frequency data and 
show no large outburst. On the other hand, the spectra before 
and after the campaign show possible outbursts, the most signif- 
icant one being the one seen by Metsahovi around MJD 55039. 
According to shock models, the outburst should move outwards 
within the jet, i.e. from high to low frequencies. As coverage 
at higher frequencies prior to the Metsahovi flare is missing, it 
cannot be excluded that a flare occurred first at these frequen- 
cies. The OVRO lower frequency observations gave a hint for a 
higher flux during their measurement four days after the flare, 
but the increase is statistically not significant due to the rather 
large uncertainty of the point. Hence no conclusions on the va- 
lidity of the shock scenario can be drawn from this data. 

The light curves measured by Effelsberg from MJD 54155 
to MJD 54952 (see Fig. [8) reveal an ambivalent correlation of 
the variability. On the one hand, there are features where only 
near-by frequency bands showed the same trends with gradu- 
ally decreasing tendency (MJDs 54239, 54486 and -54935). 
On the other hand, structures may have been present at all mea- 
sured bands with about the same strength (MJDs ~ 54547 and 
~ 54779). The different features are possibly attributed to re- 
acceleration of particles within the jet. The occurrence of an 
equal amplitude at all frequencies or gradually decreasing ampli- 
tude with frequency can be explained in that context by different 
physical conditions within the jet, e.g. a change of the magnetic 
field or the particle density. Alternatively, the sparse sampling 
in combination with frequency-dependent time lags may explain 
some of the observed features. 

OVRO and Effelsberg observe at about the same frequency 
(15.0GHz and 14.60GHz, respectively), giving the possibility 
of an independent cross-check. Indeed, the corresponding light 
curves show the same trends, and the measured data points in 
Fig. lA.ll match well. A detailed correlation analysis between the 
two instruments is not meaningful though due to the different 
light curve sampling. The well-sampled OVRO data can be used 
to overcome missing observations at 14.60 GHz around the ra- 
dio high states at MJDs 54547 and 54779 and, despite having 
some gaps at that time, indicate indeed an increased flux den- 
sity with an amplitude comparable to the one of the remaining 
frequencies. 

The probability that the flux density seen by OVRO was 
constant during this second Effelsberg high state, between MJD 
54761 and 54796, is rather low Ctf 2 /d.o.f. = 21.2/12, i.e. 
~4.8%), due to the first measurement in this time period (see 
also Fig. [5]). Neglecting this, the ^ 2 /d.o.f. = 2.8/11, giving 
highly significant evidence for constancy. The flux density rose 
from the first point within two days by ~ 16% and remained 
constant afterwards. This indicates that the peak seen in the 
Effelsberg light curve around MJD 54779 was indeed a broad 
high flux density plateau. From the OVRO variability time scale, 
the Doppl er factor can be estimated to be > 3.4 using Eqs. (1) 
and (2) in lLahteenmaki & Valtaojal (1 19991) . It should be noted 
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that this estimation method has not been tested for faint radio 
sources like TeV BL Lacs, and that the estimation of the flare 
rise time is based on two data points only. Therefore, the value 
is not representative, but also not in disagreement with the quasi- 
simultaneous results from the high state SED modelling (see 
Sect.03). 

The parsec-scale radio spectrum is fiat (see Fig. [6} show- 
ing no clear signs of the synchrotron self-absorption turnover 
at low frequencies. In principle, such a spectral shape may be 
explained as optically thin synchrotron emission from an ensem- 
ble of electrons haying a v ery hard energy spectrum N(E) oc E~ x 
dSokolovskv et al.l2010bl) . However, the more likely explanation 
is that the flat s pectrum is a result of optical ly thick synchrotron 
emission from a Bl andford & Konigl (fl979) type jet. This expla- 
nation is su pported by the ob served core size increase at lower 
frequencies dUnwin et al.l 19941 see Table[2| and the difference in 
separ ation between the comp onent C 3 and th e core observed a t 
15.4 dPiner & Edwardsll20TJ4T) and 43.2 GHz dPiner et alj|2010h . 
Together, these points agree with the standard interpretation of 
the parsec-sca le radio core in 1ES 2344 +5 14 as a surface in 
a continuous Blan dford & Konigll d 19791) jet a t which the op- 
tical depth at a given frequency v is r v =* 1 dLobanovlll998t 
ISokolovskvetal.ll2011l) . This is a challenge to most o f the alter- 
native interpr etations of the core physics discussed bv lMarscherl 
(20061, 120081) . at least for the frequency range covered by our 
VLBA observations. 

That no superluminal motion is observed in the jet of 
1ES 2344+514 is in line with the pr evious findings that this 
source and a nu mber of other TeV dPiner & Edwardsl [2004; 
iPineret alJbOld) and non-TeV dKarouzos et alJl2012h BL Lacs 
show much slower apparent jet speeds compared to those typi- 
cally found in compact extra galactic rad i o sour ces dLister et alj 
2009b). Note however, that iPiner et alJ d2010h report the de- 
tection of significant component motion in 1ES 2344+514 with 
speeds inconsistent with the results presented above. The pos- 
sible sources of this discrepancy inc lude (i) the s maller number 
of observational epochs available to IPiner et al. | ( 120101) and (ii) 
the fact that the authors combine component positions measured 
at different frequencies without explicitly taki ng the effect of 
a frequency-dependent core shift into ac count dLobanovl 1998; 
Sokol ovskv etal]|20Tll:lHada et alj|20ll) . which may introduce 
systematic errors in the positions of the components. 

There is a possibility that the observed jet component mo- 
tion is not indicative of the actual jet flow speed in this source. 
However, that assumption is not supported by the fact that the 
core brightnes s temperature is well below the inverse Compton 
limit ~ 10 12 K dKellermann & Paulinv-Tothll969HKovalev et al] 
120051) . The rather low brightness temperature of the core indi- 
cates that the radio emitting plasma in the jet is probably affected 
by only moderate relativistic beaming. 



4.2. Cross-Band Correlations and Variability Studies 

Correlated variability at different energy bands or lack of such 
provides important information on the emission mechanisms 
and locations. During the MW campaign, significant variabil- 
ity was only present in the low frequency radio and X-ray 
regime. Due to the unfortunate sampling and overall faintness 
of 1ES 2344+514, only a limited number of simultaneous inter- 
band data pairs is available, rendering it difficult to provide 
meaningful correlation studies. The most promising data sets 
collected during the campaign are the OVRO and Swift XRT 
data. During the small and short flux density rise at the be- 
ginning of the OVRO measurements the XRT flux was declin- 

































\ 






















\ 














- 


























\. _ 
























i 
























, , , 


, , , 


, , , 


, , , 


i 


0. 


4 0. 


5 0. 


6 0.17 0. 


8 



Fig. 13. Integral flux measured by Swift XRT from 2-10keV 
versus the flux density measured by OVRO for data pairs with 
a time difference of < 0.9 days. A constant fit is represented by 
the dashed line, the solid line shows the result of a linear fit. 



ing from the flare peak down to ~15%. There may, therefore, 
be an anti-correlation between these two frequency bands. In 
Fig- El the data pairs of the two instruments with a time differ- 
ence < 0.9 days are plotted against each other. Mathematically, 
a linear fit has a probability of 45.2% Or 2 /d.o.f. = 4.7/5) and 
is clearly preferred above the constant fit (logarithmic likelihood 
ratio test probability in favour of the linear fit is 98.9%). The 
hint of a correlation is dominated by the lowest flux density point 
though, and the statistical basis is weak. Therefore, no conclu- 
sive picture can be drawn from the OVRO -XRT comparison. 

Since the MAGIC light curve points are only marginally sig- 
nificant, the feasibility of investigating correlations with other 
bands is limited. It should however be noted that the daily flux 
point with the highest significance in the MAGIC light curve 
appeared only ~ two days after the highest Swift XRT flux had 
been detected. Since this high X-ray flux was accompanied by 
a strong hardening of the spectral index, this can be interpreted 
as an injection of fresh electrons into the emission region, which 
should cause a higher flux also at gamma-ray energies. 

Considering time scales beyond the duration of this MW 
campaign, flux (density) changes in the radio, optical and X- 
ray regime were clearly detected for 1ES 2344+514, as can be 
seen from the first four light curves in Fig. [14] A fit with a con- 
stant results in ^- 2 /d.o.f. = 3929.6/317 for OVRO, ^ 2 /d.o.f. = 
2058.3/437 for KVA and^ 2 /d.o.f. = 491.8/66 for Swift XRT. 
Also the distribution histograms of flux (density) over error 
(Fig. IA.31 > show a clear deviation from a Gaussian function, 
where the latter would be expected for uniformly sampled light 
curves dominated by statistical fluctuations. The strong flare 
measured by XRT around MJD 54442.2 cannot be unambigu- 
ously identified in the Effelsberg or KVA light curves. In the 
KVA light curve, a slightly higher flux was seen ~5.3 days 
after the large flare and ~5.7 days after a smaller XRT flare 
(MJD 54466.5), suggestive of a time lag of the optical emission 
with respect to the X-rays. Also the Effelsberg measurements 
revealed two significant peaks (at different frequencies, though) 
in that time period (~ MJDs 54486 and 54547). But on a signifi- 
cant correlation of these high states with the XRT flares can only 
be speculated due to the incomplete sampling of the Swift XRT, 
KVA and Effelsberg light curves. 

On time scales of years, the good sampling allows us to per- 
form a search for correlations between the OVRO and KVA data. 
In order to exclude a bias of the result caused by measurement 
noise, OVRO data with an error > 0.02 Jy (~ 15 % of the data) 
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Fig. 14. Long term light curves of 1ES 2344+514 measured by Effelsberg, OVRO, KVA+Tuorla, Swift XRT, RXTE ASM, Swift BAT and INTEGRAL ISGRI. The first four 
panels show an observation-wise binning (< 1 day), the last three panels have been re-binned to weekly scales. The blue points mark the MJD range shown in Fig. [5] The lines 
in the top panel do not result from a fit but simply connect the data points. Horizontal short-dashed lines give fits with a constant. The black solid vertical lines denote the 
arrival time of LAT events > 100 GeV. The ASM light curve has been cut at 01/01/2011. For clearness, two points of the BAT light curve are not shown (MJD 54177.5, rate 
(2.26 ± 0.91) ■ 10~ 3 cnts s _1 per detector; MJD 54541.5, rate (3.16 + 1.82) • 10~ 3 cnts s _1 per detector) and t he (linear) error bar of one point (MJD 53729.5) has been cut at the 
plot edge. Distribution histograms of the light curve points over their respective errors are given in Fig. IA.3I See text for details. 
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were excluded from this analysi s. Using the discrete cor rela- 
tion function (DCF) as defined in Edelson & Krolikl (Il988l) . we 
searched for possible correlations for lags up to +100 days be- 
tween both data sets. Two such searches were performed, one 
in which the raw light curves were used, and one in which we 
searched for correlations after first subtracting off a low-pass fil- 
tered version of the data in order to remove long-term trends 
which might influence the calculation of the DCF. The analysis 
did not yield a significant correlation. 

Investigating the light curves of RXTE ASMEl Swift BAT 
and INTEGRAL ISGRTjon time scales of one day, some outliers 
become apparent. These are expected from a statistical point of 
view, and all data points but one do not have a signal/error sig- 
nificantly offset from their corresponding Gaussian distribution 
(see Fig. IA.4I ). This point, having a signal to noise ratio of ~ 5, 
was measured by ASM at MJD 54468.0, 1.5 days after an in- 
creased flux seen by Swift XRT and ~ 3.8 days before the higher 
KVA state (see above). If real, it indicates that XRT detected the 
onset of a flare potentially even higher than the large one around 
MJD 54442.2, which preceded a small flare in the R-band by 
~ 3 - 4 days. The sparse sampling does not allow to draw further 
conclusions on the nature or origin of the flare. 

A fit with a constant to the daily light curves is ruled 
out on high statistical basis for ASM and BAT (^ 2 /d.o.f. = 
4215.6/3007 and 2364.1/1733, respectively), though not for 
ISGRI Ot 2 /d.o.f. = 228.8/214). The Gaussian fits to the sig- 
nal/error distributions reveal a significant shift of the mean 
value to positive values for BAT and ISGRI, but not for ASM 
(0.18 + 0.03, 0.11 + 0.07 and - (0.07 + 0.02) for BAT, ISGRI 
and ASM, respectively). Note, however, that a Gaussian statis- 
tical behaviour is not expected for ASM due to coded mask ob- 
servations. In the case of BAT, Gaussian statistics is still applica- 
ble despite applying the coded mask technique due to the large 
number of individual detector elements. Consequently, the BAT 
light curve indicates significant variability of the source at hard 
X-rays. 

The large flare detected by XRT on MJD 54442.2 is not 
clearly visible in the daily ASM or BAT light curve; ISGRI did 
not observe at that time. 1ES 2344+514 seems to be too faint at 
X-rays to be detected on daily time scales by these two instru- 
ments. Therefore, the light curves of ASM, BAT and ISGRI have 
been re-binned in the same way as the simultaneous BAT data 
(see Sect. 12.41 ). As an example, the weekly results are shown in 
Figs.fT4land lA.3l From the signal/error distribution histograms, 
no significant flares outside the Gaussian distributior(3 are ap- 
parent for any binning and instrument except for the 1 -week BAT 
point at MJD 54772.5 already discussed in Sect. 13.31 The large 
XRT flare at MJD 54442.2 is still not clearly present in any of 
the light curves for any binning. 



12 RXTE ASM data were obtained from NASA GSFC's archive. In 
the generation of the light curve only single-dwell ASM data were used 
in which \he x 2 red was < 1.3. Slight variations in the signal to noise ra- 
tio over the full ASM light curve are due to episodes where the source 
position was less well covered by the individual cameras of the ASM. 
Due to a strongly reduced instrument sensitivity resulting from degra- 
dation of the detectors towards the end of the mission, no data since 
01/01/201 1 have been used. 

13 Data taken from HEAVENS jWalter et alJ I201CL 
http://www.isdc.unige.ch/heavens/). 

14 A closer look at the light curves reveals several extended peaks with 
rather symmetrical rise and fall times in the ASM and BAT data for 
different binnings, but these structures match quite well the minima of 
the solar angle to ASM and can therefore most probably not be ascribed 
to 1ES 2344+514. 



The trend of a positive Gaussian mean increases with 
larger time bins for BAT, finally leading to the detection of 
1ES 2344+514 as reported in the 58-Month Catalog. The light 
curves of ASM and BAT are not consistent with a constant flux 
up to quarterly binning, though the corresponding probabilities 
are rising with increasing time bin size. For a yearly binning, a 
fit with a constant yields ^- 2 /d.o.f. = 14.9/7, ^- 2 /d.o.f. = 5.1/5 
and ^ 2 /d.o.f. = 8.2/5 for ASM, BAT and ISGRI, respectively. 

The arrival times of the nine photons with energy > 100 GeV 
detected by Fermi-LNT within the first 44 months are also shown 
in Fig. [14] No exceptional behaviour is visible in the light curves 
of the other energy bands at these times. 

At the time of the Metsahovi flare (MJD 55039, see 
Sect. 14.1.41 ). optical R-band monitoring data are also available. 
1ES 2344+5 14 showed only an insignificant hint of a higher flux 
density at MJDs 55040 and 55041. The source was detected ev- 
ery day in the R-band at low flux levels without significant out- 
bursts from seven days before the Metsahovi flare until three 
days after. Enlarging the time window to several weeks prior 
and posterior to the radio flare also does not reveal a sign of an 
R-band flare, though it should be noted that the object was no 
longer covered on daily time scales during that time window. 
Either there was no correlation present between the R-band and 
37 GHz during this flare (as the missing long-term OVRO - KVA 
correlation is also suggesting), or the optical monitoring missed 
it. A missing correlation would hint on different flaring mecha- 
nisms or, more likely, spatially separated emission regions. In the 
latter case, a time delay between the radio and optical emission 
would be expected, which may be more firmly determined based 
on continuous monitoring in the future. Quite interestingly, two 
of the Fermi-LXT events with energies > 100 GeV were de- 
tected 46 days before and < 2 days after the Metsahovi flare, 
respectively (see Table [3). On a correlation can only be specu- 
lated though, taking into account that the exact time of the flare 
maximum cannot be determined neither from the Metsahovi nor 
the OVRO or KVA monitoring data. 

In general, the present monitoring programs at various wave- 
lengths represent a major progress towards the understanding of 
blazar phenomena. Nevertheless, more efforts seem necessary, 
increasing the sampling density and time basis, and especially 
extending the monitored energy range to the X-ray and VHE 
regime. 

4.3. Spectral Energy Distribution Modelling 
4.3.1. Simultaneity 

Since the gamma-ray observations were only marginally signif- 
icant, (quasi-)simultaneous data sets for constructing SEDs are 
composed according to the X-ray flux state. We define a "low" 
and "high" X-ray flux SED, choosing MJDs 54760.9 (high) and 
54768.8 (low). The exact observation times of the different in- 
struments around these data sets are given in Table [4] 

The two data sets are too close in time to derive individ- 
ual gamma-ray results. The corresponding ULs of AGILE-GRID 
and Fermi-LAH as well as the MAGIC spectrum were aver- 
aged over the entire respective observation periods and used for 
the modelling of both SEDs. Note that no significant variability 
could be found in any gamma-ray band (see Sects. I3.1I3.2I ). the 
detection of which would exclude averaging the measurements. 

The 66-month BAT spectrum can be regarded as a mea- 
sure of the average low flux of the source, since no significant 
flares are present in the BAT light curve from daily to yearly 
scales apart from the artificially high weekly flux point at MJD 
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Table 4. Observing intervals of the SED data sets. 



Data Set 


Instrument 


Observation Time 
[MJD] 




Effelsberg 


54778.947-54778.950 




OVRO 


54769.077 


1FS 2344+514 low 


Swift 
AGILE 


54768.806-54768.948 
54770.500-54800.500" 




Fermi-LAT 


54759.941-54800.897° 




MAGIC 


54759.941-54800.897 




Effelsberg 


54756.960-54756.970 




OVRO 


54761.095 




VLBA 


54761.96 -54762.42 


1ES 2344+5 14 high 


KVA 

Swift 


54761.718 

54760.899-54760.983 




AGILE 


54770.500-54800.500° 




Fermi-LAT 


54759.941-54800.897° 




MAGIC 


54759.941-54800.897 



Notes. (fl) No detection; 95 % c.l. ULs calculated. 



54772.5. Taking into account that the low and high state XRT 
spectra have been measured in a slightly lower and higher flux 
state than the long-term average (judging from Fig.fT4l). respec- 
tively, the BAT spectrum may be considered as being quasi- 
simultaneous to these spectra. 

The variability time scale during the observations in the radio 
regime is hard to assess and differs from band to band, but in 
general large changes in flux are not expected on time scales of 
~two weeks. Therefore, some radio measurements have been 
included on a quasi-simultaneous basis. 

For the chosen low and high state dates of 1ES 2344+514, 
there were no simultaneous measurements by KVA. For the high 
state, the result from the following night is used. Taking into 
account that the KVA measurements do not show a hint of vari- 
ability throughout more than three months of observations, this 
procedure seems justified. On the other hand it should be noted 
that the simultaneous optical - X-ray data pairs are very few and 
especially missing for high X-ray fluxes, hence a higher optical 
flux during the XRT flare cannot be excluded. 

4.3.2. Model Description 

Two different leptonic SSC emission models have been applied 
to the s uch defined quasi-simu l taneou s SEDs. The one-zone 
model by Mar aschi & Tavecchiol (l2003) describes the SED com- 
pletely by nine parameters: the radius R, Doppler factor 5 and 
magnetic field B of the emission region, which contains an elec- 
tron distribution following a broken power law with index n\ for 
Train < y < Tbreak and index n 2 for y bre ak < 7 < Tmax with density 
K at Lorentz factor y = 1. y,„,„ has been fixed to values of 1 
and 4000, which represent the extreme cases of the lowest and 
a very high realisation, visualising a large part of the reasonable 
parameter range. 

In the second model (Weidinger& S panierll2010h . an elec- 
tron distribution with density K at y m ; n is being accelerated in a 
zone with radius R icc . The electrons are finally escaping the ac- 
celeration region and enter a second region, i.e. the emission re- 
gion, with radius R em where no further acceleration takes place. 
The magnetic field B and Doppler factor 5 are the same for the 
two regions. The resulting electron distribution as well as the 
spectral indices are derived self-consistently from the accelera- 
tion and cooling processes and are not determined a priori. 



4.3.3. Results 

The 1ES 2344+514 SEDs compiled from this campaign are 
shown in Fig. [15] The given simultaneous KVA and UVOT data 
have been host-galaxy corrected where necessary as well as de- 
reddened. Due to the strong host galaxy and the large uncertainty 
in its flux, the V-, B- and U-band fluxes cannot be determined 
with sufficient significance and hence are not shown. Since the 
optical data given at ASDC are not host-galaxy corrected, they 
have been omitted. The simultaneous XRT data were corrected 
for galactic absorption. MAGIC and Whi pple data were EBL 
de-absorbed using the model fromlKneiske & Dold (1 2010). the 
archival VERITAS data by the lFranceschini et al.l d2008h model. 
The AGILE UL denotes the flux > lOOMeV, the Fermi-LAT 
1FGL ULs are given between 0.1 and 0.3 as well as 0.3 and 
lGeV Note that the Metsahovi UL has a c.l. of 4<x and the 
IRAM ones 3 cr. 

A comparison of our results to archival SED data reveals that 
the source has been measured in one of the lowest flux states ever 
obtained from X-ray to VHE gamma rays. At optical and radio 
frequencies, the fluxes were at a modest level. 

For comparison, we also include published Fermi-LAT 
1FGL points on a quasi-simultaneous basis, taking into account 
that no significant variability is present also in 2FGL. Though 
all points are consistent with a simple power law distribution, 
there seems to have been a small jump between the LAT points 
at 6 and 60 GeV which is hard to describe with the applied mod- 
els. Either the 6 GeV point is rather low in flux or the 60 GeV 
point comparably high. The latter point connects rather smoothly 
to the MAGIC spectrum from 2005/6 and the VERITAS points 
from 2008, whereas putting more weight on the 6 GeV point 
the MAGIC spectrum from this campaign seems to match bet- 
ter the non-contemporaneous LAT data. Since the highest en- 
ergy point is subject to a rather large statistical uncertainty (< 10 
and < 20 events in 1FGL and 2FGL, respectively) and hence 
prone to potential short-term flux variations, our modelling is 
focusing on the 0.6-6 GeV points. In this context, the inconsis- 
te ncy of the Fermi-LA T spectrum with the VHE points reported 
m Abd oltaD d2009b) is no longer evident using the MAGIC 
points from this campaign, indicating that the VHE results de- 
rived here are more representative of the average flux state of 
1ES 2344+514. 

The 66-month BAT spectrum is an adequate extrapolation of 
the XRT high state from this campaign. The low state spectrum, 
on the other hand, would require an increase of flux with rising 
energy, which cannot be described with the SSC models applied 
here. Therefore, the BAT spectrum has been considered for the 
fit to the high state data set only. 

Note that SSC models are in general not suited to describe 
the low frequency (i.e. radio) emission. Photons of these ener- 
gies are self-absorbed in the radiation field of the SSC emis- 
sion region. The observed flux in the radio regime is probably 
produced by cooled electrons from an outer region of the jet 
which are unimportant for the mo delling of the higher frequen- 
cies dMaraschi & Tavecchiol2003h . 

Taking this into account, both models are in reasonable 
agreement with the simultaneous data of 1ES 2344+514. The 
quasi-simultaneous 1FGL points disfavour the one-zone model 
fits having a y m j n = 1 . Clearly the fits described by y m \ n = 4000 
are preferred, or in general values being closer to 4000. However, 
also for the one-zone fits with high y m ; n , a softer spectral index in 
the HE regime would improve the compatibility with the 1FGL 
data. 
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Fig. 15. Simultaneous SEDs of 1ES 2344+5 14 as derived from this campaign (black, blue and red markers) together with monitoring 
results from F-GAMMA and the Tuorla Blazar Monitoring Program (green points) as well as archival data (shown in grey). The 
black MAGIC data points represent the overall flux during the campaign, "low" and "high" denote the simultaneous data sets as 
given in Table [4] For the one-zone model, the number given in the legend represent s the value of Ymin- See text for further details . 
Archival data have been taken from the ASDC SED Builder ( |http://tools.asdc.asi.it j), |Acciari etall d20 1 lbh . lSchroedter et all d2005l) 
(Whipple) and Alb ert et all ((2007b) (MAGIC). The Planck 2 a ULs were taken from lGiommi et al] (120121) . 



The derived model parameters (see Table [5j are in good 
agreement with typical values for HBLs. Only y m j n seems 
rather high for the one-zone model, as does y max for the low 
state two-zone model. Considering the exceptional faintness of 
1ES 2344+514 across several energy bands, this concordance 
is not necessarily expected. Consequently, either the low flux 
state detected in this campaign does (still) not represent the 
"quiescent" state of the source, or the quiescent state model 
parameters do not differ considerably from the already known 
ones. The ULs on size and magnetic field strength in the dom- 
inating radio emitting reg ion derived from VLBA observations 
( Sok olovskv et al.l20 1 ()a) do not contradict the parameters of the 
blazar emission zone as given in Table [5] 

Comparing the high and low state in terms of the one-zone 
model, the latter is explained by a softer electron spectral index 
as well as a lower ybreak and y max - In the case of the two-zone 
model, the magnetic field drops consistently accompanied by 
a higher y m ; n , ybreak and y max - The parameter changes of each 
model are best explained by a change in the acceleration proper- 
ties of the non-thermal electrons, i.e. the efficiency of the under- 
lying Fermi processes drops or rises respectively (we recall that 
Tbreak is not computed self-consistently in case of the one-zone 
model). This behaviour may be caused by the e mitting volume 
leaving a standing feature along the jet (see e.g. iMarscher et all 



2008) or, more likely, due to the observation of two independent 
blobs. 

It is interesting to note that the luminosity of the inverse 
Compton component of the low (i.e. low X-ray flux) state SED 
exceeds the one of the high state SED for all applied mod- 
els. Specifically in the case of the one-zone model, this makes 
the bolometric luminosities Lboi of the two flux states com- 
parable. For y min = 1 (4000), L boUow = 10 44 7 (10 44 6 ) and 
^boi.Mgh = 10 44 8 (10 447 ) erg s _1 , which is basically identical con- 
sidering the uncertainties involved. 

Due to the differences between the two model approaches 
at sub-optical frequencies and in the hard X-ray to soft gamma- 
ray band, it is in principle possible to distinguish between the 
validity of the models. The first frequency band is covered by 
Planck, though no detection of the source has been reported i n 
The Early Release Compact Source Catalogue (lAde etalJ2011l) . 
containing the results of the first ten months of operation. Swift 
BAT, INTEGRAL IBIS as well as AGILE-GRID and Fermi-LAT 
cover the second window, but are not sensitive enough to detect 
the source during low flux states on short time scales. To ex- 
clude one of these models, a more sensitive instrumentation than 
currently available is needed. 

Within this campaign, a small shift between the synchrotron 
peak in the high and low state may be present (see Fig. H3V A 
peak estimate from the data has been obtained by fitting the opti- 
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Table 5. Model parameters. See text or reference for an explanation of the models. The used data sets and their simultaneity are 
discussed in the text. 



Reference 


This Campaign 




Albert et al. (2007b) 


Tavecchio et al. (2010) 


Acciari et al. (201 lb) 


Model Reference 


(1) 


(2) 






(3) 


(1) 




(3) 


Flux Level 


low high 


low 


high 


low 


high* 




low 


high 


B [G] 


0.07 


0.05 


0.09 


0.10 


0.08 


0.1 


0.09 


0.03 





20 


26 




8 


15 


LJ 


13 


20 


flem [10 15 Cm] 


3 4 


9 


5 




10 


4 




10 


R m [10 13 cm] 




8 


13 












K a [10 5 cirr J ] 


4.5 1.9 


0.2 


0.1 


-0.5 


-0.4 


0.3 




-0.4 


e\ 


2.3 


2.5 


2.3 




2.2 


2 


2.5 


2.3 


e 2 


3.4 3.2 


3.5 


3.3 




3.2 


3.2 




3.2 




1 or 4000 


1800 


550 


-2500 


-1500 


1 (8000) c 




-200 


Ttaeak [10 4 ] 


5 8 


10 


3 




- 15 


1 


-40 


-50 


Tmax [10 6 ] 


0.7 1.5 


6.3 


1.5 


-0.8 


-1.6 


0.7 




-2.0 



Notes. (<I) Note that K is given for the acceleration region at y rain in Weidinger & Spanier (2010), whereas the value is defined for the emission 
region at y m j n = 1 for the other two models. (6) See footnote [T5] (c) Fit has been performed on optical data not corrected for the host galaxy. See 
also footnote [T6l 

References. (D lMaraschi & Tavecchid d2003h : (2) IWeidinger & Spanierl d2010h : (3) lKrawczvnski et al] d2004l) . 



cal and X-ray SED points by a parabolic power-law in apex form 
(see e.g. Tra macere et al.l[2007r> : 



VFy 



fo • lO^ Oos.oO/vpcaO) 2 ergcirT 2 s" 1 



(2) 



where v pea k is the frequency of the synchrotron peak. The re- 
sult is shown in Fig. [16] From the fit to the low state data 
(j( 2 /d.o.f. = 6.5/8), the peak energy is determined marginally 
significant as (0.27 + 0.24) keV. In the high state case, the good- 
ness of fit is higher (x 2 /d.o.f. = 5.6/10), but the fit yields 
an insignificant peak energy of (2.8 + 4.5)keV. Additionally, a 
parabolic power law is not clearly preferred over a simple power 
law in both cases (logarithmic likelihood ratio of 79.5 % and 
95.7% for the low and high state, respectively). Consequently 
the data are insufficient to directly constrain changes of the peak 
position. As determined from the one-zone y m \ n = 4000 (two- 
zone) modelling, the peak shifted from ~0.15keV to ~ 1.7keV 
(~ 0. 1 3 keV to - 0.46 keV) by ~ 1.1 (~ 0.4) orders of magnitude 
between the low and high flux state. All these values are far from 
the extreme blazar characteristics 1ES 2344+5 14 has shown dur- 
ing high X-ray flux states. This is expected due to the rather 
small flux difference between the two states observed here. 



4.4. Comparison with Archival Campaigns 

We compare the model parameters obtained from this campaign 
with the ones from three additional broad MW data sets reported 
for 1ES 2344+514. The parameters are listed in Table|5] 



4.4.1. Archival Campaign Description 

In Albert et al. ( 2007bl). a homogeneous one-zone SSC model 
(Krawczvnski et al. 2004) was used to model a low as well as 
a high state of the source. The low state data consisted of si- 
multaneous measurements of MAGIC, KVA and an ASM UL 
taken between 08/2005 and 01/2006, where archival data by 
BeppoSAX were added as a low state X-ray spectral template. 
Data by BeppoSAX and a simultaneous UL by Whipple were 
combined with an archival Whipple high state spectrum to de- 
scribe the high stat^3 Consequently the simultaneous data set 



15 We n ote that the Whippl e high state fluxes had been adapted 
wrongly in Alber t et alJ d2007bh . 
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Fig. 16. Simultaneous KVA, Swift UVOT and Swift XRT data 
from the low (red downward triangles) and high (blue upward 
triangles) state. The solid lines describe a fit with a log-parabolic 
power law as stated in Eq.|2] R-band and UV data have been de- 
reddened, the former additionally host galaxy corrected. X-ray 
data have been corrected for galactic absorption. 



was not very constraining and the SED models rather specula- 
tiv e, as noted by the autho rs. 

iTavecchio et al.1 d2010l) were using the results from the first 
three months of ferm/-LAT observations and combined them 
with archival, non-contemporaneous data of 1ES 2344+514 
from the radio up to the VHE banc(3 They used the one-zone 
SSC model also applied in this work. 

A more rece nt cam paign, also using the model of 
iKrawczvnski et all (|2004), combined measurements by 
VERITAS, RXTE PC A, Swift XRT and Swift UVOT from 
10/2007 until 01/2008 dAcciariet al.ll201 lbl) . The low state SED 



16 We note that the optical data used for the SED modelling had not 
been host-galaxy corrected. Increasing y min from 1 to - 8000 can com- 
pensate for the missing correction. 
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was modelled using the time-averaged VERITAS measurements 
(excluding a large flare) and a representative moderate X-ray 
flux 1 -day spectrum by XRT and PCA, similar to the procedure 
used in this campaign for the low state SED. Note however 
that the MAGIC light curve from this campaign did not show 
significant variability, opposite to the VERITAS measurements 
even after the exclusion of the flare (F V ar = (34+ 16)%). 
The high state SED data set was built from the highest fluxes 
measured by VERITAS and Swift XRT together with the 
corresponding UVOT data. Though the gamma-ray and X-ray 
flares seem to have been correlated, they were separated by 
> 24 hours. Note that the flux doubling time scale of the VHE 
and X-ray flare was also ~ 24 hours, thus the true simultaneous 
fluxes could have been different by a factor of ~ 2. The model 
predicted significantly higher, in the high state by ~ one order of 
magnitude, fluxes than the 1FGL points. Since these were taken 
after the MW campaign, higher fluxes in the LAT range during 
the observations indeed cannot be excluded, despite the rather 
constant emission of the source for two years since the launch 
of Fermi. 



4.4.2. Model Parameter Comparison 

While the model parameters are in general difficult to compare 
because they are not independent of each other, we can more 
easily invest i gate ge neral trends. Note tha t the m odels shown in 
lAlbert et al.l (12007bl) and iTavecchio etalj d2010h were not con- 
strained by simultaneous measurements as well as partly af- 
fected by wrongly adapted data, thus the results have to be taken 
with care. All models indicate a rather weak magnetic field and 
high Dop pler factor, similar to values typically found for HBLs 
(see e.g. ITavecchio et al]|2010h . Variability is not explained in 
a unique way. For instance, changes in B and 5 can occur, but 
do not have to. Variability mostly arises by changes in the elec- 
tron distributions, i.e. different spectral indices or energy distri- 
butions. Particularly in the case of this campaign, the spectral 
indices are harder in the high flux state. It is interesting to note 
that the size of the emission region is often not changing between 
different flux states. That is not expected if the two states corre- 
spond to the emission o f one moving and expanding blob of elec- 
trons within the jet (e.g. lAtovan & Afiaroniarj|1999tlSi kora et al. 
1200 lb at two different times. In the case of Ulbert et alj l2007b) 
this can be understood taking into account that the low and high 
state are not c ausally connected but artificially constructed. For 
lAcciari et al.l d201 lbl) . another flare occurred between the high 
and low state SEDs, destroying causal connection. Moreover, the 
time differences between the two flux states are, also in the case 
of this campaign, too large to allow both to be caused by the 
same emitting region, at least if the emission region is visualised 
as a travelling and expanding (with ~ c) blob (which would ex- 
pand by several orders of magnitudes within one week). To al- 
low for a casual connection between the high and low state, the 
emission region has to be externally confined in some way (e.g. 
by magnetic fields). An alternative and more likely explanation 
attributes the emission to a standing shock. In this case, the size 
of the emission region would not change considerably, and vari- 
ability wou ld be caused by structu ral changes due to variations 
of the flow dTagliaferri et al.ll2008h . 

In general, the parameters found in this campaign are in good 
agreement with archival values. It should however be noted that 
the range for some parameters is rather broad, due to the inter- 
dependencies of the parameters and non-unique solutions of the 
model fits. More data and more MW campaigns in different flux 



states are needed to fully constrain the models and reduce the 
allowed parameter ranges. 

5. Summary and Conclusions 

In this paper, the results from the first MW campaign on 
1ES 2344+514 from the radio to the TeV band have been pre- 
sented and discussed, also taking into account multi-band long- 
term data. The MW observations took place from 10/2008 until 
01/2009, where the ~ 40 day long core campaign was conducted 
in October and November. The source was found at low to mod- 
est flux states at radio and optical frequencies, whereas at X- 
rays and gamma rays, the flux level was amongst the lowest ever 
reported for 1ES 2344+514. Due to this faintness, HE observa- 
tions did not result in a detection during the campaign, and the 
time-averaged VHE detection was only marginally significant. 
Nevertheless we were able to obtain a reliable VHE spectrum 
due to the large observation time of ~ 20 hours, the good event 
statistics and the source being a well-known VHE emitter. 

The VHE analysis suggested a rather hard spectral index, 
which, if real, would be opposite to the "harder when brighter" 
trend found in general for HBLs. The flux was consistent with 
constant during the campaign. At X-rays, a moderate flare was 
detected with ~ halving of the flux within several days. During 
the flare, hints for a counter-clockwise behaviour in the hardness 
ratio - integral flux plane were found, indicating that the flare 
was caused by a shock front characterised by comparable cool- 
ing and acceleration time scales. This finding was strengthened 
by the constant spectral index during the flare. Taking into ac- 
count all X-ray observations during the campaign, the spectral 
index still did not show significant variability, though a potential 
correlation between the X-ray index and the X-ray integral flux 
was visible. The evolution of the hardness ratio with the integral 
flux corroborated that "harder when brighter" trend, confirming 
findings reported in the literature for 1ES 2344+514. No signif- 
icant variability could be found at optical and UV frequencies. 
From VLBA observations, the size of the radio core could be de- 
termined or constrained at several frequencies, yielding values 
of the order of 10 17 cm. This is more than one order of magni- 
tude above the size determined from SED modelling, indicating 
a different origin of the radio and SSC emission. 

1ES 2344+514 exhibited significant variability only at low 
frequency radio and X-ray bands during the campaign. Due to 
that as well as unfortunate sampling and technical problems, the 
basis for cross-band correlations for the time of the MW cam- 
paign is too small for a meaningful investigation. For a flare ob- 
served at the end of the core campaign, indications were found 
suggesting it may have been caused by injection of fresh elec- 
trons into the jet. 

On time scales longer than this campaign, significant vari- 
ability was evident for the radio, optical and X-ray regimes 
whereas the high energy gamma-ray light curve from 2FGL was 
consistent with being constant. In contrast to the low and con- 
stant emission found by ferm/-LAT, the events with energies 
above 100 GeV detected from 1ES 2344+514 indicate that the 
source may have a comparably high flaring duty cycle. 

Different feature characteristics were found in the Effelsberg 
light curve at low radio frequencies, indicating a possible re- 
acceleration of particles within the jet. The difference between 
these features may be explained by changes in the environment 
of the particles. The behaviour of the combined long-term ra- 
dio spectra of the source gave rise to interpreting the emis- 
sion as a two-component system composed of quiescent diffuse 
emission overlaid by frequent outbursts. The signature of such 
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shocks should be traceable from higher to lower radio frequen- 
cies. One flaring event at 37 GHz was visible on the investigated 
time scales, but the expected signatures could not be found in 
the other radio bands. The OVRO measurements had some gaps 
during these days but did not show evidence of a significant flux 
increase. Two of the nine photons with an energy > 100 GeV 
were detected by Fermi-LAI around the time of that flare. A 
counterpart of this event was not found in the optical R-band 
despite rather good optical coverage, which would hint at a dif- 
ferent emission region of the 37 GHz and R-band emission if not 
due to sampling effects. The flare on its own represents a rare 
event for HBLs concerning its amplitude and time scale. A long- 
term correlation analysis between the 15.0 GHz and R-band was 
conducted, yielding no significant correlation between the two 
bands. 

The observed flat parsec-scale radio spectrum together with 
the frequen cy-dependent core si z e can be interpreted as a sig- 
nature of a Blandf ord & Konigl |l979) type jet. Analysing all 
MOJAVE observations conducted until today, no significant mo- 
tion of the three identified jet components on time scales of 
eleven years could be found, opposite to claims in previous 
publications. The apparent jet speeds of the components were 
/?app < 0.13, with the most constraining value having been 
-0.01 ± 0.02. 

Monitoring at soft and hard X-rays revealed only one sig- 
nificant individual detection, though a general trend of positive 
flux for BAT and ISGRI was apparent, leading to the detection 
of the source by BAT from 58 months of data. The individual de- 
tection, found by ASM, was coincident within a few days with a 
higher state seen by XRT and a hint for an R-band flare seen by 
KVA. Also the Effelsberg measurements showed increased ac- 
tivity around this time period. However the sampling was insuf- 
ficient for a meaningful investigation of the origin of the flare. 
The BAT light curve was significantly variable. The long-term 
trend measured by ASM did not show a hint of a positive signal. 

From the observations, (quasi-)simultaneous SEDs for a low 
and high X-ray flux state were constructed and modelled using a 
one-zone SSC as well as a self-consistent two-zone SSC model. 
Both could describe the data well, however quasi-simultaneous 
HE data posed some challenges for the modelling. In particular, 
these disfavoured the one-zone models having a y m ; n of unity, be- 
ing in general better described by the upper part of the tested pa- 
rameter range. The one- and two-zone models suggested a shift 
of the first SED peak by ~ 1 . 1 and ~ 0.4 orders of magnitude, 
respectively. Direct fitting of the combined optical and X-ray 
data did not result in a firm determination of the peak energies. 
The individual parameters retrieved from the one- and two-zone 
modelling were mostly in agreement between these two differ- 
ent model approaches for each of the two flux states. They were 
consistent with values found in archival campaigns as well as 
standard parameter ranges for HBLs. This concordance is not 
self-evident in the context of "quiescent" state emission, where 
the quiescent spectrum should be dominated by a low and con- 
stant flux component which possibly has different spectral char- 
acteristics. Either the "quiescent" state was not detected within 
this MW campaign, or the corresponding model parameters do 
not differ significantly from the typical values. The two applied 
models showed significant differences at high radio frequencies 
and in the hard X-ray to HE bands. In the future, instruments 
more sensitive in these regimes could probe the validity of the 
models. 
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Table A.l. 5w/ff XRT results. 
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Notes. (n) Swift observation ID. (A) Swift XRT exposure. <c) Integral flux between 2 and lOkeV determined by a simple power law fit from 0.3- lOkeV. (< " Spectral index determined by a simple 
power law fit from 0.3 - 10 keV. <c) Integral flux determined by a log parabolic power law fit from 0.3 - 10 keV. ^ Spectral indices determined from a log-parabola fit from 0.3 -10 keV. (s) Probability 
that the log-parabolic power law fit is preferred over the simple power law fit by means of a logarithmic likelihood ratio test. {h) Hardness ratio, defined here as counts (2 -iokeV)/counts (0 .2-ikeV)- 
(,) Observation time too short for extracting results. 
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Table A.2. Analysed MAGIC data sets and results using detec- 
tion cuts to determine the significance and open cuts for the in- 
tegral fluxes. 
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Notes. <n) If dates are given, they correspond to the day following 
the observation night. < * ) Arithmetic average of observation duration. 
<c> Effective observation tim e. (<,) Significance of the signal calculated 
according to Li & Ma ( 1983) Eq. 17. (c) Measured integral flux. ULs are 
given with 95 % c.L. We recall that the fluxes and signal significances 
were determined using different cuts. 
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J. Aleksic et al.: The Simultaneous Low State SED of 1ES 2344+514 from Radio to VHE 
Table A.3. Swift UVOT results. 



Obs. ID" 


MJD starl 


Exp." 
[ks] 


Fw 

[mag] 




F B 

[mag] 




F u 

[mag] 








■fuvwi 
[mag] 






F\]VM2 

[mag] 




^uvw 
[mag] 




35031019 


54730.158 


0.89 


15.45 


+ 0.10 


16.47 


± 0.10 


16.69 


± 





12 


17.27 


± 


0.20 


17.76 


± 


30 


17.76 


± 0.15 


35031021 


54745.554 


1.63 


15.46 


±0.10 


16.49 


±0.10 


16.57 


± 





12 


17.20 


± 


0.20 


17.79 


± 


30 


17.63 


±0.20 


35031022 


54749.513 


0.88 


15.47 


±0.07 


16.45 


±0.07 


16.53 


± 





10 


17.15 


± 


0.10 


17.57 


± 


15 


17.41 


±0.10 


35031023 


54757.769 


1.14 


15.46 


±0.10 


16.52 


± 0.10 


16.71 


± 





12 


17.30 


± 


0.20 


17.74 


± 


30 


17.49 


±0.15 


35031024 


54759.895 


2.18 


15.47 


±0.07 


16.51 


±0.07 


16.66 


± 





10 


17.31 


± 





10 


17.73 


± 


15 


17.66 


±0.10 


35031025 


54760.899 


2.26 


15.46 


±0.07 


16.52 


±0.07 


16.69 


± 





10 


17.27 


± 





10 


17.66 


± 


15 


17.62 


± 0.10 


35031026 


54761.904 


2.23 


15.46 


±0.07 


16.47 


±0.07 


16.56 


± 





10 


17.16 


± 





10 


17.68 


± 


15 


17.66 


± 0.10 


35031027 


54762.908 


2.38 


15.45 


±0.07 


16.46 


±0.07 


16.52 


± 





10 


17.17 


± 





10 


17.41 


± 


15 


17.54 


± 0.10 


35031028 


54763.167 


4.83 


15.46 


±0.05 


16.51 


± 0.05 


16.57 


± 





10 


17.19 


± 





10 


17.59 


± 


15 


17.57 


± 0.10 


35031029 


54764.857 


1.52 


15.50 


±0.07 


16.46 


±0.10 


16.50 


± 





10 


17.10 


± 





15 


17.33 


± 


20 


17.50 


± 0.10 


35031030 


54765.917 


2.48 


15.41 


±0.07 


16.50 


±0.07 


16.57 


± 





10 


17.31 


+ 





10 


17.58 


±0 


15 


17.62 


±0.10 


35031031 


54766.865 


1.15 


15.41 


±0.10 


16.49 


± 0.10 


16.61 


± 





10 


17.20 


± 





10 


17.56 


±0 


30 


17.55 


±0.15 


35031032 


54767.869 


2.68 


15.48 


±0.07 


16.51 


±0.10 


16.65 


± 





10 


17.37 


± 





10 


17.80 


±0 


20 


17.78 


±0.10 


35031033 


54768.806 


1.25 


15.44 


±0.07 


16.50 


±0.07 


16.60 


± 





10 


17.23 


± 





10 


17.66 


± 


20 


17.67 


± 0.10 


35031034 


54769.932 


1.82 


15.50 


±0.07 


16.52 


±0.07 


16.67 


± 





10 


17.48 


± 





10 


17.62 


± 


15 


17.70 


±0.10 


35031035 


54770.881 


2.05 




































35031036 


54771.933 


1.63 


15.49 


±0.07 


16.51 


± 0.07 


16.73 


± 





10 


17.38 


± 





10 


17.66 


± 


15 


17.58 


± 0.10 


35031037 


54772.892 


0.18 






16.54 


±0.07 


16.77 


± 





10 


17.22 


± 





10 












35031038 


54773.892 


0.91 


15.40 


±0.07 


16.44 


±0.07 


16.54 


± 





10 


17.12 


± 





10 


17.51 


±0 


20 


17.49 


±0.15 


35031039 


54777.483 


1.54 


15.39 


±0.07 


16.48 


±0.07 


16.54 


± 





10 


17.21 


± 





10 


17.62 


±0 


15 


17.59 


±0.10 


35031040 


54784.592 


1.03 


15.50 


±0.07 


16.52 


±0.07 


16.69 


± 





10 


17.36 


± 





10 


17.52 ± 


15 


17.68 


±0.10 



Notes. (fl) Swift observation ID. {b} Swift total exposure of all UVOT filters. 



Table A.5. Calculated HE luminosities and number of events 
above 100 GeV detected by Fenni-LAT from five HBLs. 



Source 


z" 


Index 6 


i60GeV C 


N d 


N/ 








[10 43 ergs" 1 ] 






Mrk421 


0.030 


1.77 + 0.01 


19.2 ± 1.1 


35 


18 


Mrk501 


0.034 


1.74 ±0.03 


6.95 + 0.75 


16 


10 


1ES 2344+514 


0.044 


1.72 ±0.08 


2.30 ± 0.56 


9 


9 


Mrk 180 


0.046 


1.74 ±0.08 


2.01 ±0.51 


1 


1 


1ES 1959+650 


0.048 


1.94 ±0.03 


6.00 + 0.80 


3 


4 



Notes. <fl) Red shift. m Simple power law spectral index measured 
by Femzi-LAT (Nolan et al. 2012). Note that for all source the simple 
power law is clearly preferred over a curved description of the spec- 
trum. w Luminosity at 60 GeV, de termined on the b asis of the 10- 
100 GeV photon counts reported in lNolan et al] d2012h . M Number of 
events above 100 GeV. {e) Number of events above 100 GeV scaled to 
the distance of 1ES 2344+514. 
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Fig.A.2. Distribution of events with energies > 100 GeV de- 
tected by Fermi-LAT from Mrk 421, Mrk 501, 1ES 2344+514, 
Mrk 180 and 1ES 1959+650. 
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Fig. A.l. Individual radio spectra of Fig. [7] 
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Table A.4. Gaussian Component Properties of 1ES 2344+514 measured with the VLBA at 15.4 GHz. 



Component 


Epoch 




(f 


S c 


Maj." 


Axial 


PA' 


log 7V 


I'app 








[mas] 


[deg] 


[mJy] 


[mas] 


Ratio 5 


[deg] 


[K] 


\pasyr'] 






1999.75 






128.9 


0.05 


1.00 




11.37 








1999.85 






136.4 
















2000.02 






126.9 
















2000.22 






134.4 
















2008.41 






107.7 


0.05 


1.00 




11.37 








2008.76 






97.9 


0.16 


0.49 


322 


10.59 






Core 


2008.81 






94.3 


0.22 


0.25 


316 


10.62 








2009.15 






78.9 


0.22 


0.21 


317 


10.62 








2009.42 






83.6 


0.05 


1.00 




11.24 








2009.51 






94.3 


0.17 


0.27 


322 


10.81 








2009.63 






81.7 


0.21 


0.30 


316 


10.50 








2009.94 






98.7 


0.17 


0.28 


317 


10.83 








2010.71 






100.6 


0.08 


1.00 




10.90 








2010.84 






116.9 


0.08 


1.00 




10.96 








1999.75 


0.543 


126.1 


3.2 


0.20 


1.00 




8.62 








1999.85 


0.688 


129.3 


3.1 
















2000.02 


0.774 


133.6 


1.8 


0.13 


1.00 




8.77 








2000.22 


0.491 


134.6 


5.2 


0.31 


1.00 




8.44 








2008.41 


0.401 


141.8 


10.1 


0.30 


1.00 




8.78 








2008.76 


0.630 


136.6 


7.7 


0.22 


1.00 




8.91 






C3 


2008.81 


0.714 


137.5 


11.0 


0.25 


1.00 




8.96 


-4± 7 


-0.01+0.02 




2009.15 


0.728 


138.3 


7.4 


0.20 


1.00 




8.97 








2009.42 


0.446 


135.5 


11.6 


0.15 


1.00 




9.41 








2009.51 


0.591 


140.4 


8.3 


0.14 


1.00 




9.36 








2009.63 


0.665 


138.8 


9.4 


0.32 


1.00 




8.68 








2009.94 


0.637 


136.8 


8.0 


0.20 


1.00 




9.03 








2010.71 


0.535 


140.2 


17.0 


0.34 


1.00 




8.89 








2010.84 


0.554 


138.4 


10.7 


0.16 


1.00 




9.34 








1999.75 


1.257 


135.4 


6.2 


0.46 


1.00 




8.19 








1999.85 


1.538 


135.9 


3.6 


0.13 


1.00 




9.06 








2000.02 


1.558 


135.3 


3.2 


0.29 


1.00 




8.29 








2000.22 


1.541 


138.6 


5.0 


0.44 


1.00 




8.13 








2008.41 


1.422 


137.7 


6.5 


0.55 


1.00 




8.05 








2008.76 


1.731 


140.5 


4.9 


0.43 


1.00 




8.14 






C2 


2008.81 


1.594 


140.5 


5.5 


0.40 


1.00 




8.25 


6± 7 


0.02 ± 0.02 




2009.15 


1.521 


141.1 


5.0 


0.40 


1.00 




8.20 








2009.42 


1.416 


138.8 


6.9 


0.56 


1.00 




8.06 








2009.51 


1.394 


139.5 


7.2 


0.55 


1.00 




8.10 








2009.63 


1.590 


141.4 


6.7 


0.73 


1.00 




7.81 








2009.94 


1.491 


141.3 


6.9 


0.50 


1.00 




8.16 








2010.71 


1.657 


142.3 


7.7 


0.66 


1.00 




7.96 








2010.84 


1.481 


141.5 


12.9 


0.75 


1.00 




8.07 








1999.75 


2.644 


146.1 


4.5 


0.76 


1.00 




7.61 








1999.85 


2.832 


141.9 


3.1 


0.73 


1.00 




7.48 








2008.76 


3.410 


142.9 


3.7 


0.75 


1.00 




7.54 






CI 


2008.81 


2.836 


141.4 


4.5 


1.16 


1.00 




7.24 


22 ±23 


0.06 ± 0.07 


2009.15 


2.986 


144.6 


5.5 


1.60 


1.00 




7.05 






2009.42 


3.007 


140.9 


3.2 


1.10 


1.00 




7.15 








2009.51 


2.553 


141.7 


2.7 


0.34 


1.00 




8.08 








2009.63* 


4.495 


146.8 


4.5 


1.80 


1.00 




6.85 








2009.94 


2.998 


142.3 


3.3 


1.44 


1.00 




6.92 







Notes. An asterisk (*) indicates a component not used in the fit. 

(fl) Distance from core. ib> Position angle with respect to the core. (c) Flux density. (d) Major axis of fitted component. {e) Axial ratio of fitted 
component. (/) Position angle of component's major axis. (g) Log brightness temperature. >h) Apparent jet speed. '"'/Japp = v app /c 
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Fig. A.3. Distribution of flux density, flux or rate divided by the error for the individual light curves shown in Fig. [14] for the last 
three panels including a fit with a Gaussian (shown in grey). See text for details. 
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Fig.A.4. Distribution of rate divided by the error for the individual daily-binned light curves of RXTE ASM, Swift BAT and 
INTEGRAL ISGRI including a fit with a Gaussian (shown in grey). See text for details. 
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